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INTRODUCTION 


African  Americans  (AAs)  who  have  prostatic  adenocarcinomas  (PrCas)  tend  to  have  poorer  outcomes,  overall,  than 
European  Americans  (EAs)  who  have  equivalent  PrCas.  The  causes  of  this  poorer  prognosis  are  not  understood.  In 
addition,  the  molecular  features  of  PrCa  that  may  vary  among  AAs  and  EAs  with  PrCa  have  not  been  identified 
adequately.  The  underlying  hypothesis  of  these  grants  is  that  molecular  characteristics  of  PrCas  vary  with  race  and  such 
molecular  differences  may  explain  the  more  aggressive  behavior  of  PrCas  of  AAs  compared  to  equivalent  PrCas  of  EAs. 
Of  note,  molecular  features  and/or  their  transduction  pathways  may  be  targets,  in  the  future,  for  molecularly  based 
therapies.  Thus,  a  major  goal  of  these  grants  is  to  characterize  potentially  aggressive  molecular  features  of  PrCas  from 
AAs  that  differ  from  molecular  features  of  equivalent  PrCas  from  EAs. 

Another  hypothesis  that  may  explain  differences  in  outcomes  of  PrCa  based  on  race,  is  that  PrCas  may  induce  local  field 
effects  that  affect  uninvolved  prostate  glands  and  that  these  field  effects  may  vary  with  race  and  the  Gleason  score  of  the 
PrCa.  Such  field  effects  have  been  described  for  the  adenomatous  polyposis  coli  (APC)  gene,  the  glutathione  S 
transferase-p  1  (GSTpl)  gene  and  RASSFIA,  a  putative  tumor  suppressor  gene.  Our  findings  that  human  epididymis 
protein-4  (HE4)  is  decreased  in  the  serum  of  patients  with  PrCa  compared  to  serum  of  patients  evaluated  for  PrCa  but 
with  negative  biopsies  may  be  caused  by  similar  field  effects  via  hypermethylation  of  the  promoter  of  WFDC2,  the  gene 
ofHE4. 

The  first  approach  is  to  identify  molecular  features  of  PrCas  that  are  different  between  PrCas  and  matching  uninvolved 
histopathologically  normal  appearing  prostate  glands  from  the  same  cases.  Subsequently,  the  molecular  markers  that  are 
different  between  PrCas  and  uninvolved  prostate  glands  in  AAs  will  be  compared  with  those  from  EAs.  These  molecular 
features  will  be  evaluated  with  a  focus  on  those  molecules  that  are  likely  to  be  associated  with  the  aggressiveness  of 
PrCas. 

Because  AAs  frequently  select  therapy  by  radiation  rather  than  radical  surgery,  less  tissue  has  been  available  to  study. 
PrCa  in  AA  patients.  In  addition,  patients  with  PrCas  with  higher  stage  and  Gleason  scores  are  sometimes  not  candidates 
for  radical  surgery,  independent  of  race.  To  improve  the  study  of  PrCa  in  these  categories  for  which  tissues  are  limiting  to 
research,  we  obtain  and  utilize  nitrocellulose  (NC)  blots  (tissue  prints)  of  prostate  biopsies.  Specifically,  NC  blots  are 
obtained  from  all  prostate  biopsies  of  a  case  (i.e.,  patent  evaluated  for  PrCa).  Total  RNA,  DNA  and  proteins  are  extracted 
from  NC  blots  of  selected  biopsies  based  on  the  presence  or  absence  of  PrCa  in  individual  biopsies.  Analyses  of  these 
extracts  are  at  various  molecular  levels,  e.g.,  epigenetically  regulated  (methylated)  promoters  of  genes  at  the  DNA  level. 
Also,  the  mRNA  and  microRNA  levels  and  protein  levels  can  be  analyzed.  For  example,  the  mRNA  extracted  from 
nitrocellulose  blots  is  analyzed  on  the  Affymetrix  platform  at  the  core  facility  of  Harvard  Medical  School  to  identify 
mRNAs  associated  with  PrCa.  This  discovery  approach  is  confirmed  by  RTQ-PCR  in  a  different  set  of  samples  with 
PrCa.  Separately,  the  methylation  of  DNA  is  analyzed  focusing  on  those  genes,  GSTpl,  APC,  RASSFIA,  and  WFDC2, 
for  which  preliminary  results  indicate  that  hypermethylation  field  effects  are  induced  by  foci  of  PrCa.  Proteins  identified 
to  be  differentially  expressed  in  patients  with  and  without  PrCa  are  analyzed  in  serum,  plasma,  and  urine  by  multiplex 
((Euminex  (BioPlex))  and  in  PrCa  versus  matched  uninvolved  prostate  glands  by  liquid  chromatography  mass 
spectrometry  and  by  immunohistochemistry.  Based  on  the  discovery  results,  results  are  confirmed  in  different  sets  of 
samples. 


BODY 

Administrative:  In  2014,  the  major  administrative  issues  were  the  renewal  of  IRBs  at  UAB  and  UCA.  Because  the 
collection  of  the  contracted  60  cases  of  PrCa  was  successfully  completed  at  UCA,  the  IRB  at  UCA  is  now  in  the  data 
analysis  only  category  and  no  further  cases  are  being  accrued. 

During  2014,  several  new  personnel  joined  our  project.  These  included  two  new  urologists,  Drs.  J.  Nix  and  Soroush  Rais- 
Bahrami,  a  new  urologic  pathologist.  Dr.  J.  Gordetsky,  a  new  postdoctoral  fellow.  Dr.  D.  Otali,  and  a  new  technician,  Ms. 
A.  Stewart.  Dr.  Nix,  Dr.  Otali  and  Ms.  Stewart  have  been  added  by  amendment  to  the  UAB  IRB  of  this  project.  We  have 
submitted  an  amendment  to  add  Drs.  Soroush  Rais-Bahrami  and  J.  Gordetsky  to  the  IRB. 

During  2014,  Dr.  J.  Nix,  a  new  urologist  at  UAB  joined  our  grant  as  a  collaborator.  Both  he  and  Dr.  Soroush  Rais- 
Bahrami  specialize  in  using  overlays  of  magnetic  resonance  imaging  (MRI)  over  real  time  ultrasound  images  (MRI-US) 
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to  identify  suspicious  (atypical)  areas  of  the  prostate  that  are  likely  to  harbor  foci  of  PrCa.  This  adds  critical  technology  to 
our  study  of  the  field  effects  induced  by  foci  of  prostate  cancer. 

In  2014,  Dr.  Grizzle  continued  to  be  a  consultant  to  the  DOD  Prostate  Cancer  Tissue  Repository.  In  addition,  he 
continued  as  a  participating  member  of  the  External  Advisory  Committee  of  the  DOD  Repository  of  Lung  Cancer.  Both 
these  activities  have  ended  or  are  in  the  tissue  distribution  phases  and  no  additional  patients  are  being  accrued.  Dr. 

Grizzle  is  still  available  for  consultation  to  these  activities. 

In  2014,  Dr.  Gaston  has  served  on  the  NCI  Special  Emphasis  Panels  for  the  review  of  grant  applications  for  Cancer 
Genetics  and  Innovative  Technologies  for  Cancer  Research.  During  this  last  year.  Dr.  Gaston  served  as  a  reviewer  for  the 
DOD  Prostate  Cancer  Research  Program  (Detection  Diagnosis  and  Prognosis  Peer  Review  Panel). 

Specific  Scientific  Progress  and  Results: 

Collection  of  Nitrocellulose  Blots  for  Analysis:  In  this  last  reporting  period,  66  study  subjects  were  enrolled  for  gene 
expression  studies  for  a  cumulative  enrollment  of  1 17  patients,  83  AA  and  34  EA.  Specifically,  nitrocellulose  blots  (i.e., 
tissue  prints)  were  obtained  from  each  of  the  diagnostic  prostate  biopsy  cores  obtained  primarily  via  ultrasound  (US) 
guidance  on  these  cases.  To  date,  prostate  biopsy  tissue  samples  from  all  study  subjects  have  been  transferred  to  Dr. 
Gaston’s  laboratory.  Based  on  12  nitrocellulose  blots  on  each  case,  this  results  in  996  nitrocellulose  blots  from  AAs  and 
408  NC  blots  from  EAs.  The  details  of  the  pathology  diagnosis  are  outlined  in  Tables  1  and  2.  Overall,  46%  of  our  study 
subjects  were  diagnosed  with  prostate  cancer,  21%  with  high  grade  prostate  cancer  (Gleason  sum  score  7  or  more)  and 
25%  with  low  grade  prostate  cancer  (Gleason  sum  6).  In  our  AA  and  EA  study  subjects,  the  %  cancer  diagnosis,  %  high 
grade  cancer  and  %  low  grade  cancer  were  very  similar. 

In  this  reporting  period,  a  new  technology  of  obtaining  biopsies  using  magnetic  resonance  imaging  fused  with  active 
ultrasound  (MRI-US)  has  been  incorporated  into  the  studies.  MRI-US  identifies  areas  suspicious  for  PrCa  based  on 
molecular  effects  identified  by  this  technology.  The  suspicious  areas  are  sampled  by  biopsy  much  more  intensively  than 
non-suspicious  areas  but  non-suspicious  areas  also  are  biopsied  using  standard  systematic  approaches.  The  biopsies 
obtained  by  MRI-US  are  nitrocellulose  (NC)  blotted  just  as  with  standard  US  directed  biopsies.  The  addition  of  this 
technology  to  our  project  is  permitted  by  obtaining  the  MRI-US  equipment  and  the  addition  of  two  faculty  members.  Dr. 
Jeffrey  Nix  and  Dr.  Soroush  Rais-Bahrami,  both  of  whom  have  extensive  experience  in  using  this  methodology,  and 
support  of  a  radiologist.  Dr.  John  V.  Thomas,  who  is  also  trained  in  this  technology.  MRI-US  biopsies  are  more  likely  to 
have  PrCa  and/or  uninvolved  prostate  glands  closely  adjacent  to  PrCa  than  systematic  biopsies  obtained  by  US.  Thus, 
such  samples  should  provide  more  information  on  the  PrCa  as  well  as  the  interaction  of  PrCa  with  closely  adjacent 
uninvolved  prostate  glands.  In  this  reporting  period,  biopsies  from  2  EA  and  I  AA  were  obtained  by  this  MRI-US 
technology. 

Most  features  of  sporadic  cancer  are  not  directly  caused  by  differences  in  the  germ  line  genomic  sequences.  Although 
different  patterns  of  single  nucleotide  polymorphisms  (SNPs)  may  increase  the  risk  for  developing  specific  cancers; 
instead  epigenetic  regulation  by  promoter  methylation  of  gene  expression  of  specific  proteins  and  the  post-transcriptional 
regulation  of  molecular  features  that  affect  the  development  and  progression  of  PrCa  are  more  important  (McNally  et  al 
2013,  Grizzle  et  al  2011).  One  focus  in  these  studies  is  on  small  and  large  non-translateable  RNAs  that,  for  example,  may 
regulate  mRNAs  either  to  increase  or  to  decrease  the  phenotypic  expression  of  specific  proteins. 

Analysis  of  Epigenetic  DNA  Hypermethylation  Field  Effect  Biomarkers 

During  this  last  year,  we  completed  an  important  study  in  which  tissue  prints  of  diagnostic  prostate  biopsies  were  used  to 
determine  if  the  field  effects  generated  by  cancer-associated  DNA  methylation  might  serve  as  a  useful  test  for  the 
presence  of  occult  prostate  cancer  in  the  tissues  adjacent  to  the  biopsy  core.  Many  studies  have  shown  that  tumor- 
associated  changes  in  DNA  methylation  can  act  as  biomarkers  for  the  presence  of  prostate  cancer.  We  have  identified 
three  DNA  hypermethylation  markers  (GSTpl,  APC  and  RASSFIA)  with  field  effects  that  extend  far  enough  around  the 
histological  boundary  of  a  prostate  cancer  focus  to  be  detected  in  a  biopsy  core  several  millimeters  away.  Importantly,  we 
found  that  both  the  magnitude  and  the  extent  of  field-effect  hypermethylation  are  sensitive  to  the  grade  of  the  adjacent 
prostate  cancer.  To  ensure  that  this  study  included  enough  African  American  subjects  to  determine  if  this  field  effect  test 
could  be  used  in  this  high  risk  population,  we  included  a  set  of  study  subjects  from  Alabama  as  well  as  from  New 
England. 
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Tissue  prints  were  prospectively  collected  from  biopsy  cores  from  study  subjects  undergoing  diagnostic  prostate  biopsies. 
DNA  from  all  12  cores  was  analyzed  for  the  DNA  hypermethylation  markers  glutathione-S-transferase  p  1  (GSTPl), 
adenomatous  polyposis  coli  (APC)  and  ras  association  domain  family  member  1  (RASSFl).  Levels  of  hypermethylation 
are  reported  relative  to  a  reference  gene,  beta  actin  (ACTB).  We  compared  three  groups  of  patients: 

Group  1  (Controls):  All  12  cores  negative  for  cancer  and  for  high  grade  PIN,  no  suspicious  histology. 

Group  2:  Biopsies  show  a  limited  amount  of  low  grade  (Gleason  3+3)  prostate  cancer. 

Group  3:  At  least  1  of  12  cores  positive  for  higher  Grade  (Gleason  7  or  greater)  prostate  cancer 

(See  Field  Effect  Study  Overview  -  Figure  1) 

We  observed  a  statistically  significant  increase  in  field  effect  hypermethylation  in  the  histologically  benign  cores  from 
control  (group  1)  to  group  2  to  group  3  for  each  of  the  three  markers  (table  3).  This  observation  suggests  that  this 
epigenetic  change  in  the  benign  tissue  surrounding  a  prostate  cancer  is  sensitive  to  the  Gleason  grade  of  the  tumor.  The 
level  of  hypermethylation  in  the  cancer  cores  also  shows  correlation  with  Gleason  grade.  Interestingly,  the  level  of 
hypermethylation  is  biopsy  cores  with  Gleason  6  prostate  cancer  shows  a  different  pattern  in  Group  2  (only  Gleason  6  in 
the  gland)  than  in  Group  3  (also  some  higher  grade  cancer  elsewhere  in  the  gland).  See  Table  4  and  Figure  2. 

This  prototype  “field  effect”  test  is  particularly  timely  as  more  patients  are  asked  to  consider  active  surveillance  rather 
than  immediate  treatment  based  on  a  biopsy  diagnosis  of  low  grade  prostate  cancer.  We  note  that  published  studies  report 
that  for  several  DNA  markers  the  extent  of  hypermethylation  is  higher  in  prostate  cancers  from  AAs  than  EAs  and  an 
initial  review  of  GSTpl,  APC  and  RASSFl  markers  in  a  small  group  of  self-identified  AA  patients  suggests  that  the 
extent  of  field  effect  hypermethylation  may  be  more  extensive  (more  cores  involved)  than  in  similar  EA  patients. 

Inasmuch  as  our  plan  of  work  for  DOD  project  PC093309  already  calls  for  the  preparation  of  DNA  samples  for  the 
analysis  of  ancestry  informative  markers,  we  plan  to  use  a  new  set  of  AA  and  EA  study  subjects  from  the  Birmingham 
Alabama  area  for  validation  studies  to  further  confirm  these  findings. 

Field  effects  of  PrCa  on  adjacent  uninvolved  prostate  glands:  Our  data  indicate  that  mRNA  for  the  adenomatous 
polyposis  coli  gene  (APC),  the  gene,  glutathione  S-transferase  pi  (GSTpl),  and  the  putative  suppressor  gene,  RASSFIA, 
are  suppressed  in  uninvolved  prostate  glands  associated  with  cases  with  PrCa.  The  extent  of  suppression  of  these  genes  in 
uninvolved  prostate  glands  is  thought  to  increase  the  Gleason  core  of  the  PrCa.  During  this  reporting  period,  we  used 
immunohistochemistry  to  evaluate  whether  the  suppression  of  genes  resulted  in  decreased  phenotypic  expression  of 
protein  in  the  uninvolved  glands  adjacent  to  PrCa. 

For  the  gene  GSTp  1 ,  the  GSTp  1  protein  was  not  detected  in  PrCa  in  1 8  of  20  cases  that  were  evaluated  as  demonstrated  in 
Figure  3. 

For  uninvolved  glands,  the  staining  of  GSTp  1  was  strongest  in  the  basal  cells,  in  which  the  most  prominent  staining  was 
nuclear  with  slightly  less  staining  of  the  cytoplasm.  There  was  less  staining  of  the  luminal  cells  in  which  most  staining 
was  weak  to  moderate  with  nuclear  and  cytoplasmic  about  equivalent,  but  with  slightly  stronger  membranous  and 
perinuclear  staining.  In  general,  the  staining  of  the  uninvolved  glands  for  GSTp  1  was  quite  variable  and  the  lower 
staining  of  uninvolved  glands  adjacent  to  PrCa  but  this  was  not  always  consistent;  however,  the  difference  in  GSTpl 
staining  between  uninvolved  glands  adjacent  to  PrCa  and  away  from  PrCa  was  much  less  than  would  be  expected  based 
upon  downregulation  of  the  gene  for  GSTp  1 . 

All  20  cases  of  PrCa  stained  weak  to  moderately  for  APC  as  is  demonstrated  in  Figure  4.  For  the  uninvolved  prostate 
glands,  there  was  not  as  clear  a  difference  in  staining  of  luminal  versus  basal  cells  as  was  noted  for  GSTpl.  The  strongest 
staining  in  both  basal  and  luminal  cells  was  perinuclear  followed  by  staining  of  the  cellular  membrane.  There  was  less 
staining  of  the  cytoplasm.  While  in  general  the  staining  for  APC  in  the  uninvolved  gland  adjacent  to  PrCa,  as  with 
GSTpl,  was  somewhat  less,  the  difference  was  small  and  conclusions  were  difficult  to  make  as  to  suppression  of  protein 
expression  by  PrCa. 
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Multiplex  (Luminex,  Bioplex)  Immunoassays:  The  Luminex  (BioPlex)  technology  of  immunoassays  permit 
concomitant  assays  of  up  to  100  individual  molecules  on  one  specimen  and  approximately  34  specimens  can  be  analyzed 
in  duplicate  in  the  same  assay.  Immunoassays  of  cytokines  and  tumor  biomarkers  can  utilize  serum,  plasma,  urine,  as 
well  as  lysates  of  cells  and  tissue.  Our  prior  results  comparing  serum  and  plasma  found  that  most  molecules  are  higher  in 
plasma  than  in  serum.  Unexpectedly,  in  most  cases,  the  cytokines  in  plasma  were  lower  in  patients  with  PrCa  than  in 
patients  suspected  of  having  PrCa  due  to  of  elevated  prostate  specific  antigen  (PSA)  in  serum  but  without  PCa  on  biopsy. 

Our  efforts  to  identify  racial  differences  in  PrCa  using  Luminex  (i.e.,  Bioplex)  technology  has  been  challenging  because 
only  patients  who  are  being  evaluated  for  PrCa  are  included  in  the  study.  There  are  no  normal  controls  because  such 
comparisons  do  not  represent  a  realistic  medical  evaluation  because  young  individuals  and  individuals  without  elevated 
PSA  are  seldom  evaluated  for  PrCa. 

Our  study  consists  of  two  groups  -  one  with  PrCa  identified  by  biopsy  and  the  control  group  who  do  not  have  PrCa  based 
on  biopsy.  Patients  identified  to  have  high  grade  PIN  and/or  with  histories  of  other  malignancies  are  excluded  from  both 
groups;  however,  there  remains  the  problem  that  patients  in  whom  PrCa  was  not  identified  on  biopsy  still  may  harbor 
PrCa  that  was  missed  on  biopsy  and  is  indolent  or  is  progressing  slowly.  In  each  study,  we  determine  the  current  status  of 
patients  without  PrCa,  if  available  from  the  medical  record.  Also,  some  patients  in  both  groups  may  have  undiagnosed 
malignancies  when  samples  were  obtained. 

In  our  initial  review  of  levels  of  cytokines  and  tumor  biomarkers  in  serum  and  plasma  from  patients  with  PrCa  and  cases 
evaluated  for  PrCa  in  which  no  PrCa  was  found  on  biopsy,  we  identified  that  human  epididymis  protein-4  (HE4)  was 
greatly  decreased  in  patients  with  PrCa  compared  to  patients  in  whom  no  PrCa  was  identified  in  standard  prostate  biopsies 
(Figure  5).  Because  PrCa  seldom  replaces  most  of  the  uninvolved  glands  of  the  prostate,  lower  blood  levels  in  patients 
with  PrCa  is  most  consistent  with  decreased  HE4  in  the  uninvolved  glands  of  patients  with  PrCa. 

Of  note,  HE4,  encoded  by  the  WFDC2  gene,  is  a  secreted  glycoprotein  whose  expression  is  restricted  in  most  normal 
tissues  but  is  usually  associated  with  expression  in  the  reproductive  tract  where  it  has  been  found  to  be  important  in  the 
maturation  of  sperm  but  also  is  likely  to  have  other  functions,  especially  in  females.  In  contrast  to  the  results  in  prostate 
cancer,  high  levels  of  expression  HE4  is  indicative  of  a  diagnosis  of  ovarian  papillary  serous  and  ovarian  endometriod 
carcinomas  (about  90%  of  patients  express  high  levels  of  HE4)  and/or  the  recurrence  of  these  cancers  following  therapy; 
HE4  also  is  increased  in  patients  with  endometrial  cancers.  These  results  have  a  very  strong  association  with  these 
cancers  so  they  have  been  accepted  by  the  FDA  as  diagnostically  useful  for  monitoring  ovarian  carcinomas. 

In  contrast  to  ovarian  and  endometrial  cancers,  our  results  demonstrated  decreased  levels  of  HE4  in  the  serum  of  patients 
with  PrCa  compared  to  patients  without  PrCa.  This  confusing  result  led  to  our  study  of  HE4  based  on  the  literature, 
immunohistochemistry  and  mass  spectrometry.  Our  interest  was  based  upon  the  hypothesis  that  HE4  might  be  from  a 
gene  hypermethylated  by  a  field  effect  induced  by  PrCa  similar  to  that  seen  in  GSTpl.  Alternatively,  HE4  might  be  only 
apparently  decreased  due  to  some  type  of  bias.  For  example,  the  lower  levels  could  be  secondary  to  the  decay  of  HE4  on 
storage  of  samples.  Similar  sample  effects  have  been  reported  for  cytokines  in  serum  samples  from  patients  with 
melanomas  stored  for  more  than  just  2  to  3  years  (Potter  et  al  2012). 

The  literature  indicates  that  in  contrast  to  serum  from  ovarian  carcinomas,  the  mRNA  of  HE4  (WFDC2)  in  tissues  of  PrCa 
is  strongly  decreased  compared  with  HE4  in  uninvolved  prostate  glands  (Kim  et  al  201 1). 

Consistently,  the  literature  indicates  that  the  WFDC2  promoter  is  frequently  methylated  in  cancer  (17/22  and  6/6 
malignant  cell  lines)  but  not  in  normal  prostate  tissue  (0/10,  or  the  normal  PrCa  cell  line)  (Kim  et  al  2011).  The 
methylation  of  WFDC2  in  uninvolved  prostate  glands  adjacent  to  PrCa  has  not  been  adequately  evaluated;  however,  it  has 
been  suggested  it  may  be  similar  to  GSTpl,  and  hence,  suppress  HE4  in  adjacent  uninvolved  glands. 

We  also  have  performed  immunohistochemistry  on  prostate  tissues  from  radical  prostatectomies  with  PrCa  and  prostate 
tissues  from  radical  cystectomy  specimens  without  PrCa.  Staining  of  PrCa  was  compared  to  staining  of  uninvolved 
matching  prostate  glands,  both  from  radical  prostatectomies  and  benign  prostatic  hyperplasia  and  normal  prostate  glands 
from  radical  cystectomy  specimens.  Immunostaining  of  a  case  with  PrCa  is  demonstrated  in  Figure  6.  This  pattern  of 
HE4  staining  is  typical  for  all  the  20  cases  of  PrCa  evaluated  by  immunohistochemistry.  The  immunostaining  indicates 
that  HE4  is  expressed  in  PrCa  in  the  perinuclear  area,  the  cytoplasm  and  the  cell  membranes.  In  the  uninvolved  prostate 
glands  in  cases  with  PrCa,  HE4  is  most  strongly  expressed  in  the  basal  cell  in  the  same  general  intracellular  pattern  as 
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described  for  PrCa.  When  prostates  removed  in  a  radical  cystectomy  operation  in  which  there  is  no  PrCa  were  identified 
were  immunostained  for  HE4,  the  corpora  amylacea  stained  moderately  to  strongly  for  HE4  (Figure  7).  This  would  be 
consistent  with  the  secretion  of  the  HE4  and  the  incorporation  of  this  glycoprotein  into  the  corpora  amylacea. 

Of  note,  the  cells  of  PrCa  stained  at  about  the  same  intensity  as  the  basal  cells  of  uninvolved  glands  (Figure  6).  We  next 
began  an  evaluation  of  the  expression  of  HE4  by  mass  spectrometry  in  macrodissected  PrCa  and  matching  uninvolved 
glands. 

Discovery  of  proteins  in  PrCa  using  mass  spectrometry:  In  2014,  we  continued  to  work  with  our  collaborator.  Dr. 
James  Mobley,  who  heads  the  proteomics  facility  of  the  UAB  Comprehensive  Cancer  Center.  This  core  contains  cutting 
edge  mass  spectrometry  instrumentation  (e.g.,  a  Hybrid  Orbitrap  Velos  Pro  mass  spectrometry  system  (Thermo)  with  an 
in  line  sampler  from  Agilent)).  Specifically,  this  latest  instrumentation  of  mass  spectrometry  will  be  used  to  survey  the 
global  proteome  extracted  from  the  clinical  prostate  specimens  (frozen,  FFPE,  and  bodily  fluids).  This  system  can 
identify  from  protein  digests  more  proteins  than  any  other  system  on  the  market  (about  2000).  Our  preliminary  studies 
have  confirmed  that  as  small  as  1  pg  of  protein  extracts  from  paraffin  embedded  tissues,  frozen  and  fresh  tissues,  as  well 
as  bodily  fluids  can  be  analyzed  by  this  system.  Extracted  proteins  from  such  samples  and  digested  and  cleaved  by 
trypsin  are  separated  using  a  20  cm  X  75  micrometer  in  side  diameter  column  that  separates  the  peptide/digested 
fragments  based  on  a  3  hour  gradient.  Data  analysis  incorporates  the  Refine  MS  from  gene  data  as  well  as  other  local  data 
analysis  systems. 

In  2014,  we  found  that  paraffin  embedded  biopsies  of  our  cases  could  not  be  utilized  for  discovery  analysis  based  on  mass 
spectrometry  without  compromising  the  diagnostic  material.  Thus,  for  our  initial  discovery  studies,  we  shifted  to  using 
macrodissected  PrCa  and  matching  uninvolved  prostate  glands  from  radical  prostatectomies.  In  our  pilot  study,  PrCa  and 
uninvolved  prostate  glands  macrodissected  samples  from  8  AA  patients  and  12  paired  macrodissected  PrCa  and 
uninvolved  prostate  gland  samples  from  FA  patients  were  extracted,  processed  and  prepared  to  run  on  the  Hybrid  Orbitrap 
Velos  Pro  system.  Our  initial  focus  was  to  be  on  the  molecules  SFASE,  EOF,  CEA,  IE-13,  G-CSF  and  MIP-lp  which 
previously  were  found  to  be  elevated  in  both  the  serum  and  plasma  from  patients  confirmed  by  biopsy  to  have  had  PrCa 
when  samples  were  compared  to  the  control  group.  Of  less  interest  are  those  molecules  previously  found  by  our  BioPlex 
studies  to  be  decreased  in  serum  and  plasma  from  PrCa  cases  compared  to  controls  including  HE4,  il-la,  il-lb,  il6,  and  il- 
8. 


Our  initial  mass  spectrometry  evaluation  of  paired  macrodissected  PrCa  and  uninvolved  prostate  glands  from  8  AA  and  12 
EA  patents  (40  samples  total),  identified  that  more  protein  was  required  for  optimal  analysis  by  mass  spectrometry  than 
could  be  extracted  from  the  unstained  sections  provided  to  the  proteomic  laboratory  in  many  of  these  cases.  Successfiil 
data  from  MS  were  obtained  from  4  EA  pairs  of  specimens  and  from  5  AA  pairs  of  specimens.  Evaluation  of  the  protein 
extraction  indicated  that  extraction  of  protein  from  paraffin  embedded  tissues  is  much  less  efficient  than  extraction  of 
protein  from  frozen  tissues  and  some  of  the  extracted  material  caused  problems  with  the  mass  spectrometry  system  (e.g., 
caused  loss  of  sensitivity  of  detectors  requiring  frequent  cleaning  of  the  MS  system).  On  the  cases  for  which  inadequate 
protein  was  extracted,  additional  paraffin  sections  have  been  provided  to  the  Proteomic  Core  and  proteins  are  now  being 
extracted  from  the  new  samples.  These  results  confirmed  our  view  that  paraffin  embedded  biopsies  cannot  be  used  in 
these  discovery  studies.  For  adequate  data  analysis  from  the  specimens,  it  is  optimal  to  analyze  the  MS  data  on  all 
specimens  together.  Thus,  we  are  waiting  until  the  8  additional  EA  cases  and  4  additional  AA  cases  are  extracted  and 
successfiilly  run  on  the  mass  spectrometry  instrumentation. 

Gene  Profiling  Studies 

Affymetrix  Whole  Transcriptome  mRNA  Gene  Expression  Profiling: 

This  year,  we  have  successfiilly  used  Affymetrix  HTA  2.0  arrays  to  analyze  gene  expression  profiles  of  tissue-print  RNA 
obtained  from  biopsies  from  AA  and  EA  men  who  were  diagnosed  with  high  volume,  high  grade  prostate  cancer.  This  is 
an  important  project  milestone  that  addresses  a  significant  disparity  in  the  representation  of  AA  patients  with  high  grade 
disease  in  prostate  cancer  gene  expression  datasets.  Patients  with  extensive  high  grade  cancer  on  biopsy  are  usually 
treated  by  hormonal  therapy  and/or  radiation,  not  surgically,  and  as  a  result  only  biopsy  tissue  is  available  for  biomarker 
studies.  Moreover,  in  the  Birmingham  area  as  in  many  other  parts  of  the  country,  AA  patients  are  more  likely  to  choose 
radiation  over  surgery  for  prostate  cancer  treatment,  which  has  further  limited  the  representation  of  AA  in  biorepository 
collections  that  depend  on  radical  prostatectomy  specimens. 
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We  have  chosen  to  use  the  Affymetrix  Human  Whole  Transcriptome  Array  2.0  as  the  most  up  to  date  and  appropriate 
expression  array  for  our  biomarker  discovery  studies.  Our  strategy  for  the  whole  transcriptome-discovery  study  is 
illustrated  in  Figure  8.  For  our  AA  and  EA  study  groups  diagnosed  with  high  grade  prostate  cancer,  we  are  focusing  on 
subjects  with  biopsy  cores  that  contain  50%  or  more  high-grade  carcinoma.  Our  study  design  compares  RNA  expression 
profiles  from  high  grade  and  high  %  cancer  cores  with  RNA  from  tumor-adjacent  cores  from  the  same  prostate  that  are 
not  involved  with  cancer  (or  with  the  least  involved  core,  if  all  cores  are  cancer-positive).  Also,  we  are  also  comparing 
RNA  expression  from  our  subjects  with  high  grade  prostate  cancer  with  RNA  from  biopsies  of  “no  cancer”  control 
subjects  all  of  whose  biopsies  show  no  cancer,  no  HGPIN,  no  atypia  and  no  notable  inflammatory  changes. 

In  our  first  set  of  Affymetrix  gene  expression  studies  we  are  utilizing  samples  from  20  study  subjects:  5  AA  and  5  EA 
with  high  grade  prostate  cancer  and  5  AA  and  5  EA  “no  cancer”  controls.  .  As  with  most  gene  profiling  techniques,  the 
HTA  2.0  array  performs  best  with  high  quality  RNA.  As  shown  in  our  QC  studies  summarized  in  Figure  9,  we  routinely 
prepare  RNA  from  biopsy  tissue  prints  with  RINs  of  over  7  (consistent  with  high  RNA  quality)  and  the  total  RNA  per 
prostate  biopsy  print  is  approximately  200  ng.  It  should  be  noted  that  we  have  sufficient  RNA  remaining  from  each  biopsy 
core  to  perform  qrtPCR  confirmatory  studies  of  the  Affymetrix  array  “hits”. 

Our  biopsy  RNA  profiling  studies  are  designed  both  to  discover  individual  biomarkers  associated  with  high  grade  prostate 
cancer  in  AA  and  EA  patients  and  to  identify  molecular  subgroups  of  prostate  cancer  that  may  be  more  prevalent  in  men 
of  African  or  European  ancestry.  It  should  be  noted  that  the  Affymetrix  HTA  2.0  array  is  well  suited  for  this  type  of 
analysis;  it  is  currently  the  most  comprehensive  array  for  interrogating  human  transcripts  for  expression  profiling  and 
includes  both  coding  and  long  non-coding  RNA  transcripts.  In  addition  to  gene-level  detection,  the  HTA  2.0  provides  the 
necessary  coverage  and  accuracy  required  to  detect  all  known  human  transcript  isoforms  produced  from  a  gene. 

We  are  in  the  early  stages  of  our  HTA  2.0  array  studies,  and  we’re  exploring  different  approaches  to  differential  gene 
expression  analysis.  The  recent  publication  by  Gorlov  et  al  (2014)  has  provided  some  particularly  helpful  guidance  for 
this  project.  This  bioinformatics  group  notes  that  although  a  typical  approach  to  analyzing  tumor  gene  expression  is  to 
compare  cancer  to  adjacent  uninvolved  tissue,  an  analysis  of  inter- individual  tumor-to-tumor  variation  in  gene  expression 
is  often  a  more  efficient  way  to  identify  genes  that  are  over  or  under  expressed  in  a  molecular  subgroup.  One  important 
advantage  to  this  type  of  analysis  is  that  it  is  not  confounded  by  tumor  associated  changes  in  adjacent  normal-looking 
tissue  (cancer  “field  effects”).  In  the  biopsies  from  some  of  our  study  subjects  with  extensive  prostate  cancer,  there  are  no 
histologically  benign  cores  and  the  approach  of  Gorlov  et  al  provides  a  useful  alternative  analytical  approach. 

An  example  comparing  the  two  approaches  to  differential  gene  expression  analysis  in  our  prostate  biopsy  samples  is 
shown  in  Figure  10.  Our  pairwise  analysis  of  tumor-to-tumor  variation  in  biopsies  from  AA  and  EA  patients  with  high 
grade  prostate  cancer  has  identified  several  robustly  overexpressed  “outlier”  genes  of  interest.  These  include  fatty  acid 
binding  protein  5  (FABP5),  which  has  previously  been  identified  as  an  “outlier”  in  gene  profiling  studies  and  defines  a 
PrCa  molecular  subgroup  with  a  prevalence  of  30-40%  (Alshalalfa  et  al  2012).  FABP5  overexpression  has  been 
implicated  in  a  more  aggressive  prostate  cancer  phenotype  (lymph  node  metastasis)  (Pang  et  al  2010).  Of  particular 
interest  is  an  association  between  FABP5  overexpression  and  enhanced  tumor  angiogenesis  in  PrCa  xenograft  mouse 
models  (Adamson  2003).  Prostate  cancer  induced  changes  in  blood  flow  and  tissue  permeability  can  be  imaged  non- 
invasively  by  MRI,  and  new  technologies  that  allow  fusion  of  MRI  images  with  the  ultrasound  used  to  guide  prostate 
biopsy  are  moving  into  the  clinic.  With  the  addition  of  Drs  Nix  and  Rais-Bahrami  to  our  study  team,  we  are  planning  pilot 
studies  to  determine  if  prostate  cancers  that  overexpress  FABP5  represent  an  aggressive  subtype  that  can  be  differentiated 
by  prostate  MRI. 

Studies  Related  Studies  to  Support  This  Grant: 

We  continue  to  work  in  a  very  successful  collaboration  with  Dr.  Clayton  Yates  at  Tuskegee  University  focusing  on  racial 
differences  in  microRNAs  in  prostate  cancer.  Our  laboratories  continue  to  work  with  the  laboratory  of  Dr.  Yates  in 
characterizing  racial  differences  in  microRNAs  of  PrCas.  From  September  2013  to  August  2014,  22  paired  samples  (22 
PrCA  and  22  matching  uninvolved  tissues)  from  8  AAs  and  12  similarly  paired  samples  of  PrCa  and  matching  uninvolved 
prostate  from  EAs  were  macrodissected  and  total  RNA  (mRNA  and  microRNA)  was  extracted  from  the  macrodissected 
tissue.  The  RNA  was  provided  to  Dr.  Yates’  laboratory  for  confirmation  of  the  prior  identification  of  miRs  which  are 
expressed  differentially  in  PrCa.  Previously,  1 1  miRs  were  found  to  be  expressed  differentially  between  PrCa  and 
uninvolved  prostate  glands  at  statistically  significant  p  values. 


Analysis  of  the  new  tissues  for  microRNAs  will  be  used  to  confirm  the  prior  results  including  the  results  that  miR- 1 52  is 
correlated  with  increased  metastasis  and  biochemical  recurrence,  probably  via  downregulation  of  DNA  methyltransferase 
(Theodore  et  al  2014). 

Future  Directions: 

The  identification  by  our  Affymetrix  studies  that  the  expression  of  fatty  acid  binding  protein  5  (FABP-5)  at  the  mRNA 
level  may  be  a  marker  for  an  aggressive  subset  of  PrCas  is  an  important  observation;  FABP-5  has  been  reported  to  be 
increased  in  cell  lines  of  PrCa  as  well  as  PrCa  tissues  and  has  been  hypothesized  to  be  a  tumor  suppressor  gene  (Das  et  al 
2001);  however,  this  and  related  observations  were  not  pursued  and  we  will  verify  this  observation  at  the  protein  level. 
Because  increased  levels  in  plasma  of  FABP-5  as  well  as  FABP-4  are  biomarkers  of  metabolic  syndrome  (Bagheri  et  al 
2010),  a  validated  ELISA  assay  is  available  for  FABP-5  (Biovendor  Laboratory  Medicine).  We  will  use  this  ELISA  assay 
to  compare  plasma  levels  of  FABP-5  from  patients  with  PrCa  versus  without  PrCa  and  in  patients  with  PrCa,  the  levels  of 
FABP-5  in  plasma  for  different  Gleason  scores  of  PrCa. 

As  part  of  the  study  of  FABP-5,  immunohistochemistry  will  be  performed  on  prostates  with  PrCa  and  prostates  without 
PrCa  (radical  cystectomy  specimens).  The  intensity  of  immunostaining  for  FABP-5  will  be  correlated  with  the  Gleason 
score  to  characterize  more  specifically  the  subsets  of  PrCa  which  expresses  FABP-5.  The  cases  with  varying  staining  will 
be  analyzed  by  RT-Q-PCR  to  verify  the  mRNA-protein  association. 

Although  we  were  initially  confused  about  our  results  on  HE4  and  other  down  regulated  proteins  in  patients  with  PrCa 
compared  to  those  without  PrCa,  we  have  now  accepted  that  low  level  mRNA  of  the  WFDC2  gene  and  the  decreased 
matching  HE4  protein  may  be  molecular  markers  of  PrCa.  Thus,  in  the  next  12  months,  we  will  verify  the  HE4  results 
using  uniquely  designed  multiplex  assays.  These  assays  will  include  biomarkers  which  our  preliminary  data  indicated  are 
both  increased  and  decreased.  Paired  samples,  PrCa  and  controls,  will  be  matched  in  analysis  based  on  sample  storage 
times  (See  Challenges  and  Problems). 

Early  in  the  next  funding  period,  the  pilot  mass  spectrometry  study  will  be  analyzed  and  its  results  will  be  completed  and 
correlated  with  the  results  from  other  levels  of  analysis  (Luminex,  Affymetrix,  RT-Q-PCR  assays).  Based  on  these 
discovery  studies,  a  list  of  molecules  of  interest  in  aggressive  prostate  cancers  will  be  developed  and  proteins  will  be 
extracted  from  the  nitrocellulose  blots  to  confirm  the  results  from  paraffin  tissues. 

We  will  continue  to  work  with  Dr.  Clayton  Yates  in  the  analysis  of  microRNAs  in  AAs  versus  EAs.  The  initial 
microRNA  differentially  expressed  in  PrCa  will  be  continued  in  additional  macrodissected  cases.  As  with  microRNA  152 
(See  References  -  Theodore  et  al  2014),  the  importance  of  each  microRNA  in  PrCa  will  be  confirmed  and  racial 
differences  will  be  evaluated.  Also,  additional  tissues  will  be  macrodissected  to  separate  PrCa  from  uninvolved  prostate 
glands  from  radical  prostatectomies  of  the  prostate.  Total  RNA  will  be  extracted  and  provided  to  Dr.  Yates’  laboratory 
and  unstained  sections  for  extraction  of  proteins  will  be  provided  to  Dr.  Mobley  for  analysis  by  mass  spectrometry.  As  in 
prior  cases,  AA  and  EA  cases  will  be  matched  to  detect  racial  differences.  In  addition,  unstained  sections  will  be  provided 
to  Dr.  Gaston’s  laboratory  for  analysis  by  RT-PCR  analysis  for  racial  admixture  studies. 

Candidate  biomarkers  that  arise  from  either  the  RNA  or  protein  discovery  platforms  are  ranked  based  on  the  strength  of 
association  with  Gleason  grade,  strength  of  association  with  racial  ancestry  and  evidence  of  likely  biological  relevance 
(either  from  published  studies  or  pathway  analysis).  We  evaluate  all  top  ranked  marker  candidates  in  an  independent 
series  of  samples  using  qrtPCR,  immunohistochemistry  or  other  appropriate  methods.  For  the  set  of  markers  that  are 
confirmed  in  an  independent  series  of  samples,  we  will  focus  our  confirmatory  studies  on  the  “hits”  that  appear  to  define 
prostate  cancer  subgroups  that  are  sufficiently  prevalent  to  be  useful  clinical  markers.  For  our  Affymetrix  HTA  2.0 
GeneChip  studies,  our  primary  focus  will  be  on  defining  RNA  expression  patterns  that  correlate  Gleason  grade  and  racial 
ancestry,  including  a  detailed  comparison  of  gene  expression  profiles  in  cancers  with  a  similar  histology  from  AA  and  EA 
patients. 

For  these  studies  we  are  working  closely  with  Dr.  Towia  Liebermann,  who  is  the  director  of  the  Harvard  Cancer  Center 
Genomics  laboratory  where  we  perform  our  Affymetix  GeneChip  analysis.  Dr.  Liebermann  is  a  well-recognized  cancer 
research  scientist  and  an  expert  in  gene  expression  analysis.  Our  expression  studies  include  a  systems  biology  analyses 
that  integrates  divergent  data  into  a  comprehensive  analysis  of  key  pathophysiological  pathways  and  analyses  that 
compare  the  prostate  cancer  GeneChip  data  generated  by  this  project  with  other  publically  available  prostate  cancer  gene 
expression  datasets  to  further  characterize  potentially  important  subgroups.  In  addition,  our  data  analysis  for  the 
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Affymetrix  HTA  2.0  array  gene  discovery  studies  includes  a  comparison  of  each  PrCa  RNA  profile  with  a  series  of 
baseline  reference  RNA  samples  from  “no  cancer”  prostate  tissues  to  help  control  for  confounding  field  effects  in  benign 
tissue  adjacent  to  cancer.  We  anticipate  that  analysis  of  gene  expression  in  “no  cancer”  prostate  tissues  from  AA  and  EA 
men  who  are  free  of  prostate  cancer  on  biopsy  may  help  identify  racial  differences  that  are  associated  with  cancer  risk. 

In  the  next  reporting  period,  we  will  complete  several  manuscripts  focusing  on  the  results  of  this  project.  This  includes  a 
manuscript  on  epigenetic  regulation  of  GSTpl,  APC  and  RASSAF  in  prostate  cancer  based  on  NC  blots,  a  manuscript  on 
HE4  (WEFC2)  and  HE4  isoforms  in  prostate  cancer,  and  a  manuscript  of  fatty  acid  binding  protein-5.  Also,  we  will 
complete  a  manuscript  on  the  use  of  NC  blots  in  the  study  of  racial  differences  in  prostate  cancer. 

KEY  RESEARCH  ACCOMPLLISHMENTS 

•  We  have  utilized  the  Affymetrix  gene  array  to  identify  a  candidate  biomarker,  fatty  acid  binding  protein-5,  which, 
based  on  the  literature,  is  associated  with  the  aggressiveness  of  tumors  in  general  and  prostate  cancer  specifically. 
The  differential  expression  of  this  biomarker  has  been  confirmed  in  a  separate  set  of  PrCas  using  RT-Q-PCR.  We 
will  confirm  these  associations  also  at  the  protein  level. 

•  We  have  concluded  the  HE4  protein  and  its  gene,  WFDC2,  are  a  potential  pair  of  biomarkers  that  are 
characteristic  of  prostate  cancer. 

•  We  have  successfully  completed  the  contract  with  Urology  Centers  of  Alabama  to  obtain  nitrocellulose  hots  on  60 
African  American  men  biopsied  for  prostate  cancer.  We  have  demonstrated  that  at  UAB  we  can  use  magnetic 
resonance  imaging  combined  with  ultrasound  imaging  (MRI/US)  to  identify  prostate  biopsy  sites  and  that 
nitrocellulose  blots  can  be  obtained  from  biopsies  during  MRI/US  procedures. 

•  We  have  successfully  macrodissected  20  radical  prostatectomy  cases  (12  EA  and  8  AA).  PrCa  was  separated 
successfully  from  uninvolved  prostate  glands.  Sections  were  provided  for  discovery  studies  using  mass 
spectrometry  (MS).  Initial  results  were  obtained  from  5  AA  tissue  pairs  (PrCa  and  uninvolved  glands)  and  4  EA 
pairs.  Initial  protein  from  the  remaining  1 1  cases  was  insufficient  for  analysis  by  MS  and  the  remaining  8  EA  and 
3  AAs  have  been  recut  and  additional  tissues  provided  for  analysis  by  MS.  Also,  microRNA  was  extracted  and 
provided  from  these  cases  to  the  laboratory  of  Dr.  Clayton  Yates  for  confirmation  and  expression  of  the 
previously  published  results  as  to  microRNA  in  prostate  cancer  (Theodore  et  al  2014).  In  addition,  tissues  from 
these  same  cases  were  provided  to  Dr.  Gaston  for  study  of  fatty  acid  binding  protein-5  and  for  analysis  of  racial 
admixtures. 

REPORTABLE  OUTCOMES 

•  The  60  cases  of  PrCa  from  AAs  have  been  collected  at  UCA  completing  their  contract  requirements.  Because  all 
required  cases  of  PrCa  have  been  collected  at  UCA,  the  IRB  for  UCA  has  entered  the  data  analysis  category  and 
no  additional  patients  are  being  accrued  at  UCA. 

•  A  preliminary  study  of  racial  differences  in  PrCa  at  the  microRNA  level  has  been  published  (Theodore  et  al 
2014).  Additional  RNA  on  matching  PrCa  and  uninvolved  prostate  from  an  additional  12  EA  and  8  AA  patients 
has  been  provided  to  Dr.  C.  Y ates  to  confirm  the  initial  observations  and  to  expand  the  analysis  of  additional 
microRNAs  identified  in  the  preliminary  study. 

•  Several  review  manuscripts  associated  with  this  grant  are  published  and/or  in  press  (listed  in  body  and  included  in 
appendix). 

•  In  2014,  all  IRBs  have  been  renewed  at  UAB  and  UCA.  TU  work  has  been  classified  as  exempt. 

CHALLENGES  AND  PROBLEMS: 

Identifying  biomarkers  of  prostate  cancer  that  are  associated  with  aggressive  subtypes  has  been  a  definite  challenge  due  to 
the  numbers  of  prostate  cancer  that  behave  indolently,  slowly  progress,  that  may  not  be  identified  on  initial  biopsy,  as  well 
as  the  observation  that  many  of  the  more  aggressive  prostate  cancers  in  EAs  tend  to  develop  in  an  older  population,  and 
hence,  may  not  be  diagnosed  until  the  tumors  are  late  stage.  Such  cases  cause  problems  in  that  cryptic  PrCas  may  be  in 
the  control  group.  Thus,  there  may  be  an  unrecognized  overlap  between  PrCa  and  control  groups.  Also,  the  samples  to  be 
analyzed  may  degenerate  on  storage  after  one  or  two  years  and  this  may  lead  to  bias  in  studies.  These  problems  are  being 
addressed  in  our  future  studies. 
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CONCLUSION 


In  this  reporting  period,  we  have  identified  using  the  Affymetrix  gene  array  of  mRNA  extracted  from  NC  blots  that  fatty 
acid  binding  protein-5  (FABP5)  is  differentially  expressed  at  the  mRNA  level  in  a  subgroup  of  prostate  cancers  and 
confirmed  this  observation  by  analysis  with  RT-Q-PCR  using  a  separate  set  of  cases  of  PrCa.  These  results  are  important 
in  that,  in  the  literature,  FABP5  has  been  associated  with  aggressive  features  of  PrCa.  The  original  descriptions  of  FABP5 
were  not  confirmed  nor  expended.  Flowever,  we  consider  that  FABP5  represents  a  biomarker  which  should  be  further 
studied  and  clinically  utilized.  Also,  the  approach  used  to  identify  and  to  verify  FABP5  as  a  biomarker  will  be  used  with 
the  other  Affymetrix  data  to  identify  additional  biomarkers  of  PrCa  aggressiveness.  Similarly,  we  consider  epigenetic 
control  of  uninvolved  prostate  tissue  by  PrCa  to  be  important  in  detecting  PrCa.  Our  observations  that  FIE4  may  be  an 
epigenetically  downregulated  molecule  will  be  explored  further.  Of  special  interest  is  differential  phenotypic  expression 
of  FIE4  in  serum.  We,  together  with  Dr.  Clayton  Yates,  will  confirm  our  reported  observations  of  microRNA  expression 
in  PrCa  and  expand  these  observations  by  determining  the  biology  in  PrCa  of  the  microRNAs  identified  to  be 
differentially  expressed  in  PrCa.  Thus,  in  this  reporting  period,  we  have  begun  to  identify  successfully  biomarkers  that 
may  be  used  clinically  to  identify  aggressive  subtypes  of  prostate  cancer  as  well  as  to  detect  cryptic  foci  of  prostate  cancer 
that  may  be  missed  based  on  routine  approaches  to  biopsy. 
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SUPPORTING  DATA 


Table  1;  Subjects  Enrolled  in  Prostate  Biopsy  Tissue  Print  Study  as  of  August  31,  2014 


All  Subjects 

All  races 

AA 

EA 

N 

%  of  group 

N 

%  of  group 

N 

%  of  group 

All  Subjects 

117 

100% 

83 

100% 

34 

100% 

No  cancer  In  any  core* 

60 

51% 

44 

53% 

16 

47% 

High  Grade  PrCa 

32 

27% 

22 

27% 

10 

29% 

Low  Grade  PrCa 

25 

21% 

17 

20% 

8 

24% 

*  May  have  HGPIN  or  other  potentially  premalignant  findings 


No  Cancer  Subjects 

All  races 

AA 

EA 

N 

%  of  group 

N 

%  of  group 

N 

%  of  group 

No  cancer  in  any  core* 

60 

100% 

44 

100% 

16 

100% 

No  Ca,  HGPIN,  suspicious 

47 

78% 

37 

84% 

10 

62% 

HGPIN,  suspicious 

13 

22% 

7 

16% 

6 

38% 

*  May  have  HGPIN  or  other  potentially  premalignant  findings 


Low  Grade  PrCa  Subjects 

All  races 

AA 

EA 

N 

%  of  group 

N 

%  of  group 

N 

%  of  group 

Low  Grade  PrCa 

25 

100% 

17 

100% 

8 

100% 

AS  Candidate** 

14 

56% 

13 

76% 

1 

12% 

Not  an  AS  candidate** 

11 

44% 

4 

24% 

7 

88% 

**UCSF  Criteria  for  AS 


High  Grade  PrCa  Subjects 

All  races 

AA 

EA 

N 

%  of  group 

N 

%  of  group 

N 

%  of  group 

High  Grade  PrCa 

32 

100% 

22 

100% 

10 

100% 

Highest  GS  3+4 

20 

62% 

16 

73% 

4 

40% 

Highest  GS  4+3  or  more 

12 

38% 

6 

27% 

6 

60% 
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Table  2 


Pre  Biopsy  Serum  PSA  Aii  Subjects 

PSA  (ng/mi) 

Aii  races 

AA 

EA 

Mean 

31.2 

41.1 

7.1 

SD 

210.0 

249.1 

7.4 

Median 

5.4 

5.4 

4.6 

SE 

19.4 

27.3 

1.3 

Lowest  PSA 

0.7 

0.7 

1.4 

Highest  PSA 

2195.0 

2195.0 

42.0 

Pre  Biopsy  Serum  PSA  minus  2  outliers 


(outliers:  2195  and  634  ng/ml) 

PSA  (ng/ml) 

All  races 

AA 

EA 

Mean 

7.2 

7.2 

■ 

SD 

6.7 

6.5 

■ 

Median 

5.4 

5.4 

4.6 

SE 

0.6 

0.7 

1.3 

Lowest  PSA 

0.7 

0.7 

1.4 

Highest  PSA 

42.0 

40.0 

42.0 
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APC  rskjs 
DatasummaiY 


Hist  0- Benign 

N 

Mean 

so 

Median 

SE 

Min 

Mac 

Group  1 

240 

14.4 

29.9 

2.5 

19 

0.0 

200.6 

Group  2 

452 

44.2 

85.4 

9.0 

4.0 

ao 

638.1 

Groups 

289 

7&8 

198.7 

22.5 

11.7 

ao 

2741.5 

GSTPi  ratios 
Data  summary 


Hist  0- Benign 

N 

Mean 

so 

Median 

SE 

Min 

Mac 

Group  1 

240 

a9 

4.1 

ao 

as 

ao 

46.1 

Group  2 

451 

3.4 

129 

ao 

ae 

ao 

165.3 

Groups 

288 

116 

55.3 

00 

12 

00 

734.3 

RASSFraios 
Data  summary 


Career  Cores 

N 

Mean 

SO 

Median 

SE 

Min 

Mac 

Group  1 

240 

410 

79.5 

06 

5.1 

0.0 

479.9 

Group  2 

A52 

900 

1829 

24.8 

86 

00 

1573.4 

Groups 

289 

124.0 

2610 

S27 

15.3 

00 

2506.4 

Dunn'S  Mu  iti  pie  Comparisorts  Test 

Comparison 

Mean  Rank 
Difference 

P  Value 

SignificatKe 

Groups  Histo  Bn 

Gioup2  Histo  Bn 

-141.8 

<0.001 

•«« 

Groups  Histo  Bn 

Groups  Histo  Bn 

-257.3 

<0.001 

••• 

Group2  Histo  Bn 

Group  3  Histo  Bn 

-11548.0 

<0.001 

••• 

Dunn's  Mu  Iti  pie  Com  parsons  Test 


Compaison 

Mean  Rank 

Difference 

P  Value 

Significance 

Groups  Histo  Bn 

Groups  Histo  Bn 

-109.2 

<0.001 

••• 

Group  1  Histo  Bn 

Groups  Histo  Bn 

-187.8 

<0.001 

••• 

Groups  Hiso  Bn 

Groups  Hiso  Bn 

-785 

<0.001 

••• 

Dunn's  Mu  Iti  pie  Com  par  isons  Test 

Comparison 

Mean  Rank 
Difference 

P  Value 

Signifkanoe 

Group  1  Histo  Bn 

Group2  Histo  Bn 

-115.4 

<0.001 

••• 

Groups  Histo  Bn 

Groups  Histo  Bn 

-148.0 

<0.001 

••• 

Group!  Histo  Bn 

Groups  Histo  Bn 

-3262 

p>ao5 

ns 

Table  3.  APC,  GSTPi  and  RASSFl  Methylation  Ratios  in  the  Histologically  Benign  Cores 


14 


APC  ratios 

Data  summary 

Dunn's  Multiple  Comparisons  Test 

Comparison 

Mean  Rank 

Difference 

PValue 

Sienificance 

Cancer  Cores 

N 

Mean 

SD 

Med  ian 

SE 

Min 

Max 

Groups  Gls 6  Group 3  Gls 6 

-28713 

<ao5 

• 

Group  2  Gis6 

&4 

114.5 

149.8 

59.8 

187 

0 

6CB.8 

Group  2  Gls  6  Group  3  Gls  3*4 

-3a722 

<aoi 

•• 

Groups  Gls6 

S7 

2214 

258.9 

1187 

42.6 

86 

1047.7 

Group  2  Gls  6  Group  3  Gls  4*3^ 

-19.117 

X15 

ns 

Groups  61s S^4 

4i 

1914 

175.7 

1386 

27.4 

3.4 

8602 

Group  3  Gls  6  Group  3  Gls  3*4 

-8009 

>a5 

ns 

Groups  6ls4»S<> 

16 

1517 

174.6 

780 

486 

17.8 

674.2 

Groups  Gls 6  Group 3  Gls 4+3*^ 

6.596 

>a5 

ns 

Groups  Gls 3*4  Group 3  Gls 4*3'* 

116 

>a5 

ns 

GSTPi  ratios 

Data  summary 

Dunn's  Multiple  Comparisons  Test 

Comparison 

Mean  Rank 

Difference 

PValue 

SiKnrficarxe 

Cancer  Cores 

N 

Mean 

SD 

Median 

SE 

Min 

Max 

Groups  Gls 6  Groups  Gls 6 

-38439 

<001 

•  • 

Groups  6is6 

64 

25.1 

50.4 

81 

6.3 

0 

283.8 

Group  2  Gls  6  Group  3  Gls  3*4 

-2848 

<006 

• 

Groups  Gls6 

37 

sao 

98.9 

310 

183 

0 

3381 

Group 2  Gls 6  Groups  Gls 4*3-* 

-28445 

>05 

ns 

Groups  GlsS*'4 

41 

8a4 

124.2 

281 

19.4 

0 

4180 

Groups  Gls 6  Groups  Gls 3*4 

4.960 

>05 

ns 

Groups  Gls4^S* 

16 

92.9 

129.2 

23.5 

32.3 

0 

437.9 

Groups  Gls 6  Group 3  Gls 4*3* 

4.994 

>05 

ns 

Groups  Gls 3*4  Group 3  Gls 4*3* 

0.034 

>05 

ns 

RASSF  ratios 

Data  summary 

Dunn's  Multiple  Comparisons  Test 

Comparison 

Mean  Rank 

DiffererKe 

PVaiue 

Si^rficarxe 

CancerCores 

N 

Mean 

SD 

Median 

SE 

Min 

Max 

Groups  Gls 6  Group 3  Gls 6 

-2175 

>05 

ns 

Groups  61s 6 

64 

27Q1 

358.4 

124.2 

44.8 

0 

1GQ82 

Group  2  Gls  6  Group  3  Gls  3*4 

-27.24 

<005 

• 

Group  S  Gls  6 

37 

419.1 

390.6 

2686 

64.2 

a4 

14480 

Groups  Gls 6  Group 3  Gls 4*3* 

-2833 

>05 

ns 

Groups  GlsS*4 

41 

515.9 

552.3 

282.7 

882 

0 

2302.6 

Groups  Gls 6  Group 3  Gls 3*4 

-850 

>05 

ns 

Groups  GIsA^S** 

16 

686.4 

983.08 

292.7 

245.77 

0 

31988 

Groups  Gls 6  Group 3  Gls 4*3* 

-4.58 

>05 

ns 

Groups  Gis3*4  Group 3  Gls 4*3* 

0.91 

>05 

ns 

Table  4:  APC,  GSTPi  and  RASSFl  Methylation  Ratios  in  the  Cancer  Cores 
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Figure  1:  Field  Effect  Study  Overview 
Primary  Question: 

Is  the  Epigenetic  Hypermethylation  of  APC,  GSTPi  and  RASSFl  Sensitive  to  the  Gleason  Grade 
of  an  Adjacent  Prostate  Cancer? 

Study  Groups: 

Group  1:  20  prostate  biopsy  patients  with  no  prostate  cancer,  HGPIN  or  suspicious  histology 

Group  2:  46  prostate  biopsy  patients  diagnosed  with  minimal  Gleason  3+3  prostate  cancer 

Group  3:  36  prostate  biopsy  patients  diagnosed  with  at  least  one  core  of  Gleason  7  (3+4  or  4+3)  prostate 

cancer 

Total:  102  study  subjects  with  12  biopsy  cores  for  each  subject 
1224  total  biopsy  cores 

Biomarker  Assay: 

Multiplex  methylation  sensitive  PCR  for  APC,  GSTPi  and  RASSFl  promoter  hypermethylation 
Major  findings: 

The  level  of  hypermethylation  in  the  histologically  benign  biopsy  cores  is  positively  correlated  with  the  level  of 
hypermethylation  in  the  cancer  cores. 

The  level  of  hypermethylation  in  the  benign  biopsy  cores  increases  significantly  with  the  Gleason  grade  in  the  cancer 
cores 

The  level  of  hypermethylation  in  a  Gleason  6  cancer  is  significantly  higher  in  Group  2  (Gleason  6  only)  than  in  Group 
3  (which  has  a  high  grade  focus  elsewhere  in  the  gland). 
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Figure  2:  Hypermethylationmarkerprofilesof  Gleason  3+3  Biopsy  Cores  in  Group  2  vs  Gleason  3+3  Biopsy  Cores  in 
Group  3 .  The  data  is  plotted  asquartiles,  with  darker  shades  of  each  color  corresponding  to  higher  levels  of 
hypermethylation.  For  all  three  markers,  the  Gleason  3+3  cores  from  group  3  are  more  highly  methylated  than  the 
Gleason  3+3  cores  from  group  2 

Group  2:  subjects  whose  cancer  cores  show  only  Gleason  3+3 
Group  3:  subjects  who  have  one  or  more  Gleason  7+  cores 


Figure  3  -  Due  to  the  large  size  of  Figure  3  and  the  email  limitations,  Figure  3  will  be  uploaded  thru  UAB’s  dropbox. 


Figure  4  -  Due  to  the  large  size  of  Figure  4  and  the  email  limitations,  Figure  4  will  be  uploaded  thru  UAB’s  dropbox. 
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Figure  5: 
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Figure  6  -  Due  to  the  large  size  of  Figure  6  and  the  email  limitations,  Figure  6  will  be  uploaded  thru  UAB’s  dropbox. 


Figure  7  -  Due  to  the  large  size  of  Figure  7  and  the  email  limitations,  Figure  7  will  be  uploaded  thru  UAB’s  dropbox. 
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Study  Groups  and  Samples 


AA 

EA 

High  Grade 

High  Grade 

Cancer 

Cancer 

AA 

EA 

No  Cancer 

No  Cancer 

Sample$:TwoRNA samples  from  each  high  grade 
cancersub)ect,onefroma  cancer  and  onefroman 
uninvolved  biopsy  core 

Cancer  Core 

•  High  Grade  Cancer  in  50%  or  more  of  the  core 

*  At  least  10,000  cells  harvestedon  the  nitrocelulose 
Uninvok/ed  Core 

*  Uninvolved  (or  least  involved)  with  cancer 

•  At  least  10,000  cells  harvestedon  the  nitraelulose 


Samples:  OneRNA sample  from  each  no<ancer 
subject,  composed  ofequalalquotsof  RNAfrom 
each  of  the  12  biopsy  cores 

No  cancer 

NoHGPIN,  no  atypia 

No  notable  inflamnutory  changes 

Pre-biopsy  PSA  nmninully  elevated 


Design  of  Affymetrix  Array  2.0  mRNA  profiling  study 

Note  that  in  some  cases  there  is  no  biopsy  core  that  is  uninvolved  with  cancer.  In  such 
instances,  we  compare  the  most  and  least  involved  cores. 


Sample  ID 

Status 

RIN  number 

RNA  cone  in 
sample  (ng/pl) 

Total  RNA  in 
sample  (ng) 

1012-A 

Tunxjr  1 

7.3 

8.06 

403 

1012-B 

Benign  1 

11 

6.11 

306 

1016-A 

Tumor  1 

7.9 

23.46 

1,173 

1016-B 

Tumor  2 

6.0 

9.95 

497 

10 14- A 

Tumor  1 

8.6 

35.14 

1,757 

1014-B 

Benign  1 

7.1 

10.11 

505 

1028-A 

Tumor  1 

~1 

28.68 

1,434 

1028-B 

Benign  1 

in 

33.99 

1,699 

2016-A 

Tunxjr  1 

3.8 

6.13 

307 

2016-B 

Benign  1 

7.0 

4.89 

244 

1007-A 

Tumor  1 

s7 

10.46 

523 

1007-B 

Tumor  2 

9.6 

19.25 

963 

1023-A 

Tumor  1 

8.3 

14.99 

750 

102  3-B 

Tumor  2 

8.3 

26.02 

1,301 

2035-C 

Entirely  Benign 

85 

9.56 

478 

202  3-C 

Errtirety  Benign 

7.6 

19.49 

975 

471-C 

Entirely  Benign 

8.7 

33.87 

1,694 

451-C 

Entirely  Benign 

7.9 

136.35 

6,817 

Figure  8  Agilent  Bioanalyzer  analysis  of  biopsy  tissue  print  RNA 
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Sample  ID 

Prosate  Biopsy  Histology 

Total  Cells  on 
Biopsy  Print 
(Quantifiler) 

RNA 

Concentration 

(ng/pl) 

Total  RNA  (ng) 

RIN 

Number 

P 1007-7 

90%  GletBon  4*5 

214,929 

19.3 

962.6 

9.6 

P  1023*4 

80%  Gleason  3*5 

161,322 

15.0 

749.6 

83 

Figure  9.  Examples  of  Bioanalyzer  Analysis  of  Biopsy  Tissue  Print  RNA  Quality 

Examples  of  biopsy  tissue  print  RNA  obtained  from  biopsy  tissue  prints  from  biopsies  with  a  high 
%  of  high  grade  prostate  cancer  are  shown  above.  Note  that  sample  P1007-7  was  prepared 
using  a  protocol  that  recovers  microRNAs  while  sample  P1023-4  was  prepared  using  our 
standard  AllPrep  RNA  DNA  protocol.  The  difference  in  recovery  of  microRNAs  can  be  visualized 
in  the  Bioanalyzer  tracing  (in  the  zone  marked  "small  RNAs") 
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Venn  Diagram 
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Affymetrix  HTA  2.0  Array 
Profiling  of  Prostate  Biopsies 
from  Patients  with  High 
Volume,  High  Grade  Cancer: 
Summary  of  global  RNA 
expression  patterns  focused 
on  differential  expression 
between  same-patient  tumor 
and  benign  biopsies 


Benign  cores  Paired  High  Grade/High  96  Cancer  Core  and  Uninvolved  (or 
No  cancer/low  risk  Least  Involved)  Core  from  the  Same  Subject 

Subjects  Park  Blue  Cancer,  Light  Blue  Uninvolved 

Expression  analysis  focused  on  Tumor-to-Tumor  comparisons  (Orlov  et  al  2014)  is  a  method  for 
identification  of  prostate  cancer  subgroups  that  minimizes  confounding  field  effects  in  tumor 
adjacent  benign  tissue.  This  has  proven  to  be  an  efficient  strategy  for  identifying  candidate 
genes  that  are  confirmed  by  qrtPCR  in  independent  samples 


Figure  10.  Tumor-to*AdJacent  Benign  and  Tumor-to-Tumor  Comparisons  of  Gene  Expression 
Profiles  for  Prostate  Cancer  Biomarker  Discovery. 
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ABSTRACT 

miRNA  expression  in  African  American  compared  to  Caucasian  PCa  patients  has 
not  been  widely  explored.  Herein,  we  probed  the  miRNA  expression  profile  of  novel 
AA  and  CA  derived  prostate  cancer  cell  lines.  We  found  a  unique  miRNA  signature 
associated  with  AA  cell  lines,  independent  of  tumor  status.  Evaluation  of  the  most 
differentially  expressed  miRNAs  showed  that  miR-132,  miR-367b,  miR-410,  and  miR- 
152  were  decreased  in  more  aggressive  cells,  and  this  was  reversed  after  treatment  of 
the  cells  with  5-aza-2'-deoxycytidine.  Sequencing  of  the  miR-152  promoter  confirmed 
that  it  was  highly  methylated.  Ectopic  expression  of  miR-152  resulted  in  decreased 
growth,  migration,  and  invasion.  Informatics  analysis  of  a  large  patient  cohort 
showed  that  decreased  miR-152  expression  correlated  with  increased  metastasis 
and  a  decrease  in  biochemical  recurrence  free  survival.  Analysis  of  39  prostate  cancer 
tissues  with  matched  controls  (20  AA  and  19  CA),  showed  that  50%  of  AA  patients 
had  statistically  significant  lower  miR-152  expression  compared  to  only  35%  of  CA 
patients.  Ectopic  expression  of  miR-152  in  LNCaP,  PC-3,  and  MDA-PCa-2b  cells  down- 
regulated  DNA  (cytosine-5}-methyltransferase  1  (DNMTl)  through  direct  binding  in 
the  DNMTl  3'UTR.  There  appeared  to  be  a  reciprocal  regulatory  relationship  of  miR- 
152/DNMTl  expression,  as  cells  treated  with  siRNA  DNMTl  caused  miR-152  to  be 
re-expressed  in  all  cell  lines.  In  summary,  these  results  demonstrate  that  epigenetic 
regulation  of  miR-152/DNMTl  may  play  an  import  ant  role  in  multipie  events  that 
contribute  to  the  aggressiveness  of  PCa  tumors,  with  an  emphasis  on  AA  PCa  patients  . 


INTRODUCTION 

The  most  commonly  diagnosed  type  of  cancer 
among  men  in  the  US  is  prostate  cancer  (PCa),  which 
accounts  for  29%  (241,740)  of  all  new  cancer  cases. 


Although  the  number  of  new  cases  of  PCa  has  decreased 
in  recent  years,  there  are  still  racial  and  ethnic  differences 
in  PCa  epidemiology.  African-Americans  (AAs)  have  the 
world’s  highest  incidence  of  PCa  and  more  than  twofold 
higher  mortality  rate  compared  with  Caucasian  Americans 
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(CAs)  [1].  Overall,  AA  patients  are  younger  and  have 
higher  Gleason  scores,  PSA  levels,  and  incidence  of 
palpable  disease  [2].  Various  factors  have  been  associated 
with  the  more  aggressive  prostate  tumors.  For  example, 
differential  gene  expression  in  AA  patients  contributes  to 
aggressive  disease  [3-7],  and  epigenetic  mechanisms,  such 
as  DNA  methylation,  result  in  the  loss  of  key  regulatory 
genes  [8,  9].  This  is  particularly  evident  for  AA  patients 
where  hypermethylation  of  genes  in  normal  or  pre- 
malignant  areas  are  thought  to  predispose  to  malignancy 
[10,  11].  However,  the  underlying  mechanism  of  these 
acquired  methylation  patterns  is  poorly  understood. 

MicroRNAs  (miRNAs)  are  small  RNA  molecules 
consisting  of  19-23  nucleotides  that  regulate  various 
biological  processes.  More  than  60%  of  protein-coding 
genes  may  be  targeted  by  miRNAs  [12],_mainly  through 
translational  repression  and  degradation  of  target 
mRNAs.  An  expanding  body  of  evidence  supports  a  role 
for  miRNAs  in  disease  progression  and  the  potential 
for  epigenetic  mechanisms,  such  as  DNA  methylation, 
to  regulate  miRNA  expression  [13,  14].  Recently,  DNA 
methylation  of  proximal  CpG  islands  in  miRNA  promoters 
was  described  as  a  method  for  decreased  expression 
in  various  cancers,  including  PCas  [15-18].  Although 
miRNAs  are  expressed  differently  in  healthy  tissues 
and  cancers  [19]  and  in  localized  and  advanced  tumors 
[20],  little  is  known  about  racial  differences  in  miRNA 
expression.  Identification  of  unique  miRNAs  and  mRNA- 
associated  targets  will  begin  to  clarify  the  specific  events 
involved  in  the  progression  of  PCas  in  AAs. 

To  address  this  question,  our  laboratory  has 
established  non-malignant  and  malignant  cell  lines 
derived  from  AA  PCas  that  replicate  many  of  the  clinical 
features  of  the  original  PCas  [21,  22].  Thus,  utilizing 
these  cell  lines  and  commonly  available  PCa  cell  lines, 
we  explored  the  possibility  of  race-related  differences  in 
miRNA  expression.  Herein,  we  report  that  AA  cell  lines 
have  a  distinct  miRNA  signature,  independent  of  tumor 
status.  Evaluation  of  the  miRNAs  most  differentially 
expressed  demonstrated  that  miR-132,  miR-367b,  miR- 
410,  and  miR-152  were  decreased  in  more  aggressive 
cells  and  that  this  expression  was  reversed  after  treatment 
of  the  cells  with  5-aza-2'-deoxycytidine  (5-aza-2’d). 
Bisulfite  conversion  and  sequencing  of  the  promoter 
showed  that  miR-152  was  highly  methylated  in  LNCaP 
and  PC-3  cells.  Ectopic  expression  of  miR-152  resulted 
in  decreased  cell  proliferation,  migration,  and  invasion.  In 
a  panel  of  39  PCa  tumors  with  adjacent  matched  controls 
(20  AAand  19  CA) ,  50  %  of  AA  patients  and  35%  ofCA 
patients  demonstrated  statistically  significant  lower  miR- 
152  expression  compared  to  adjacent  controls.  Ectopic 
expression  of  miR-152  down-regulated  DNA(cytosine- 
5)-methyltransferase  1  (DNMTl)  through  direct  binding 
in  the  DNMTl  3’UTR.  This  appeared  to  be  an  inverse 
relationship,  as  cells  treated  with  DNMTl  siRNA  re¬ 
expressed  miR-152.  Finally,  as  determined  with  a  large 


patient  cohort,  loss  of  miR-152  was  related  to  poor  clinical 
outcomes  including  decreased  biochemical  recurrent  free 
survival. 

RESULTS 


miRNA  Profile  of  Panel  of  AA  and  CA  Cell  Lines 

To  determine  the  miRNA  expression  pattern  in  our 
panel  of  AA  and  CA  prostate  cell  lines,  the  expression  of 
662  miRNAs  was  analyzed  utilizing  Asuragen  Affymetrix 
Gene  Chips.  To  validate  the  malignant  and  non- 
malignant  status  of  this  cell  line  panel,  we  first  preformed 
hierarchical  clustering  (Figure  la)  and  a  PLS  plot  (Figure 
lb),  using  reported  non-malignant  and  malignant  status  as 
the  grouping  criteria.  Both  analysis  demonstrated  distinct 
expression  patterns  in  the  malignant  and  non-malignant 
cell  lines,  with  21  differentially  expressed  miRNAs 
between  the  groups  (Supplemental  Table  1  (p<0. 00001). 
To  further  determine  race  related  miRNAs,  we  performed 
a  separate  hierarchical  clustering  analysis  on  the  same  data 
set  utilizing  racial  background  as  the  grouping  criteria. 
Distinct  hierarchical  clustering  of  miRNAs  was  evident 
between  the  AA  and  CA  cell  lines  (Figure  Ic),  regardless 
of  malignancy  status.  A  PLS  plot  reflected  this  pattern, 
with  distinct  variations  in  the  racial  profiles  of  the  AA  and 
CA  cell  lines  (Figure  Id).  In  this  variation,  there  were  47 
miRNAs  differentially  expressed  by  race  (Supplemental 
Table  2).  To  confirm  these  47  miRNAs  are  race  related, 
we  interrogated  a  microarray  dataset  that  consisted  of 
primary  tumor  cells  isolated  from  5  AA  and  4  CA  patients 
(Supplemental  Table  3).  These  race  related  miRNA  probes 
were  able  to  differentiate  patients  relevant  to  race  using 
PLS  plot  test  (Supplemental  Figure  1). 

To  if  these  miRNA’s  have  a  relationship  to 
prostate  cancer  progression,  took  the  most  significantly 
differentially  expressed  race  related  miRNA’s  (p<0. 00001) 
(Table  1),  and  performed  qRT-PCR  validation  within  our 
prostate  cancer  cell  line  progression  model.  Of  the  11 
miRNAs,  we  confirmed  the  5  most  significant  by  qRT- 
PCR  in  non-malignant  and  selected  aggressive  malignant 
cell  lines  (Supplementary  Figure  2).  With  the  exception  of 
miR-363,  which  did  not  exhibit  a  clear  expression  pattern, 
miR-132,  miR-376b,  miR-410  and  miR-152  exhibited  a 
decreased  expression  profile  as  the  metastatic  capability 
of  cell  lines  increased.  Thus,  we  chose  these  miRNAs  for 
further  analysis. 

Demethylation  Treatment  Reverses  miRNA 
Expression  of  AA- Associated  miRNAs. 

Previous  reports  and  in  silica  analyses  have 
demonstrated  that  a  majority  of  the  miRNAs  we  found 
associated  with  race  contain  CpG  islands  within  the 
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Table  1:  Hierarchical  clustering  of  race  related 
miRNAs 

miRNA 

(Adjusted)  p-value  ] 

hsa-miR-363 

3.43e-05 

hsa-miR-132 

0.000354 

hsa-miR-376b 

0.00188 

hsa-miR-410 

0.00291 

hsa-miR-152 

0.00471 

hsa-miR-189 

0.00498 

hsa-miRllO 

0.00514 

hsa-miR-27b 

0.00657 

hsa-miR-519c 

0.0071 

hsa-miR-520h 

0.00774 

hsa-miR-27a 

0.00837 

**miRNAs  separated  by  race  are  listed  based  on  p-value 

promoter  regions  upstream  of  the  start  site  (Supplemental 
Figure  3).  To  determine  if  hypermethylation  is  associated 
with  decreased  expression  of  these  miRNAs,  LNCaP  and 
PC-3  cells  were  treated  with  5  pM  5-aza-2’d  alone  for 
three  days  or  in  combination  with  100  nM  TSA  for  24  hr. 
Re-expression  of  multiple  miRNAs  was  evident  in  both 
LNCaP  and  PC-3  cells  after  5-aza-2’d  treatment.  However, 
miR-376b  and  miR-152  showed  the  most  substantial 
increases  in  both  cell  lines  (Figure  2). 

Since  miR-152  contained  the  greatest  percentage 
of  methylated  CG  sequences  and  demonstrated  the  most 
consistent  increases  after  5-aza-2’d  treatment,  we  focused 
on  this  miRNA.  To  determine  the  methylation  status  of 
the  miR-152  promoter  region,  we  extracted  DNA  from 
LNCaP  and  PC-3  cell  lines  and  preformed  sodium  bisulfite 
modification  prior  to  sequencing.  To  confirm  these 
results,  sodium  bisulfite-converted  DNA  was  subjected 
to  capillary  electrophoresis  and  analyzed  utilizing  the 
online  Bisulfite  Sequencing  DNA  Methylation  Analysis 
(BISMA)  sequencing  program.  The  results  indicated 
that  DNA  of  both  LNCaP  and  PC-3  cell  lines  was  100% 
methylated  at  1000  base  pairs  upstream  from  the  promoter 
region  (Figure  3).  Thus,  these  results  suggest  that  in 
malignant  PCa  cell  lines,  miR-152  is  inactivated  through 
hypermethylation. 


PCA  Pkt  |«7«) 


Fig.  1:  Heat  Map  of  miRNA  microarray,  (a)  Hierarchical  clustering  analyses  of  miRNAs  (left)  grouped  by  non-malignant  and  malignant 
status:  non-malignant  cell  lines  (red  block),  malignant  cell  lines  (blue  block),  (b)  PLS  plot  of  AA  and  CA  prostate  cell  lines  grouped  by 
malignancy,  (c)  Hierarchical  clustering  analyses  of  miRNAs  (left)  grouped  by  race.  AA  cell  lines  (red  block),  CA  cell  lines  (blue  block). 
miRNAs  are  ordered  according  to  their  cluster  detennined  by  p-values  using  the  Kruskal- Walls  test,  (d)  PLS  plot  of  AA  and  CA  prostate 
cell  lines  grouped  by  race. 
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miR-152  Expression  Correlates  with  Clinical  and 
Pathological  Variables. 


b. 


Fig.  2:  Multiple  mlRNA  expression  pattern  after  treatment 
with  demethylation  agent  (a)  LNCaP  cells  were  treated  with  5 
|iM  5-aza-2’d  for  4  days  alone  (+)  or  with  100  nM  TSA  for  24 
hr  (++).  Results  shown  are  representative  of  three  independent 
experiments  ±  s.e.*P<0.05.  (b)  PC-3  cells  were  treated  with  5  pM 
5-aza-2’d  for  4  days  alone  (-I-)  or  with  100  nM  TSA  for  24  hours 
(++).  Results  shown  are  representative  of  three  independent 
experiments  ±  s.e.  *P<0.05. 


The  expression  of  miR-152  in  28  normal  cell  lines, 
97  primary  tumors,  and  13  metastases  was  analyzed  using 
the  Taylor  et  al.  GSE21032  data  set  available  on  the  GEO 
website  (http://www.ncbi.nlm.nih.gov/geo/).  Primary  and 
metastatic  tumors  had  lower  levels  of  miR-152  relative 
to  normal  samples  (Figure  4a),  which  correlated  with  the 
higher  incidence  of  metastatic  samples,  metastatic  events, 
and  lymph  node  invasion  (Figures  4b  and  c).  Tumors 
with  low  miR-152  levels  also  had  reduced  biochemical 
recurrence-free  survival  (Figure  4c).  Together,  this 
describes  a  consistent  picture  of  low  miR-152  levels 
associated  with  PCa  metastasis  and  recurrence. 

These  results  prompted  us  to  determine  miR- 
152  expression  in  our  patient  cohort  of  AAs  and  CAs, 
hypothesizing  that  tumors  from  AAs  would  have  lower 
miR-152  than  those  from  CA  patients.  Patients  were 
selected  based  on  cancers  of  higher  total  Gleason  score 
(>6)  and/or  pathological  stage  (pT2),  and  positive 
for  perineural  and/or  vascular  invasion,  as  these 
pathological  characteristics  correlate  positively  with 
tumor  aggressiveness  and  metastasis.  AA  patients  had  a 
lower  median  age  relative  to  similarly  staged  CA  patients, 
which  was  also  associated  with  lower  miR-152  levels 
measured  by  qRT-PCR  (Table  2,  p<.001).  Analysis  of 
miR-152  expression  in  individual  tumors  compared  to 
matching  adjacent  nonnal  controls  showed  a  statistically 
significant  decrease  in  miR-152  expression  in  50%  of  the 
AA  patients  compared  to  only  35  %  of  CA  patients  (Figure 
5a,  b  Supplemental  table  4). 


niaaisiaziMniin 
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I  2  3  4  5  i  7  e  3  14  11  12  13  14  IS  16  17  13  13  24  21  22  23  24  25  26  27  23 
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p(2-3  X  R  R 
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Fig.  3:  Bisulfite  Sequencing  of  CpG  islands  in  the  miR-152  promoter,  (a)  Schematic  map  of  the  CpG  islands  in  the  miR-152  promoter 
using  UCSC  Genome  Browser.  CpG  islands  are  indicated  as  vertical  lines  on  map.  (b)  Bisulfite  modification  was  performed  on  LNCaP 
and  PC-3  cells,  and  DNA  was  sequenced  to  detennine  the  methylation  status  of  the  promoter  region.  Bisulfite  sequencing  was  analyzed  by 
BDPC  software  (see  Methods).  Open  square;  unmethylated  CpG;  closed  circle,  methylated  CpG. 
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Table  2:  Clinical  Characteristics  of  Paired  Primary  Tumor  Prostate  Samples 

Ethnicity 

African  American 

Caucasian 

Pt 

Samples 

20(51.3%) 

19  (48.7%) 

Age 

Median 

57.5 

64 

Mean 

57.3 

60 

Standard  deviation 

6.68 

7.3 

Min  -  Max 

47-69 

48-70 

Initial  Clinical  Stage 

T1 

NA 

NA 

T2 

16  (80%) 

10(52.6%) 

T3 

4  (20%) 

8  (42.1%) 

NS 

NA 

1  (5.3  %) 

Gleason  Score 

6 

8  (40%) 

8  (42.1%) 

7 

12  (60%) 

10  (52.6%) 

8 

NA 

1  (5.26) 

miR-152  Expression 

Average  All  Samples 
±SD 

L6± 

2.6± 

p<.00D 

tP  value  with  is  deemed  significant  utilizing  E-  test . 

NA-  none 

SD-  Standard  Deviation 

Restoring  miR-152  Expression  Decreases  Cell 
Growth. 

The  clinical  relevance  of  miR-152  in  PCa  prompted 
us  to  determine  if  loss  of  expression  had  a  biological  or 
functional  role  in  promoting  metastatic  tumors.  After 
optimizing  the  concentration  of  miR-152  mimics  that 
restored  miR-152,  miR-152  was  transfected  into  LNCaP, 
PC-3,  and  MDA-PCa-2b  cells.  As  determined  by  MTT 
assays,  ectopic  expression  of  miR-152  inhibited  cell 
proliferation  after  3  days,  and  this  continued  to  Day  6 
(Figure  6  a,b,c).  Since  the  inhibition  of  cell  proliferation 
began  at  72  hr,  miR-152  transfected  cells  were  next 
assayed  after  72  hr  by  flow  cytometry  to  determine  if 
there  was  an  effect  on  cell  cycle  progression.  miR-152 
treatment  caused  cells  to  accumulate  at  the  G2-M  phase 
(LNCaP  NC  22.90%  compared  to  miR-152-transfected 
cells  30.2%;  PC-3  NC  19.9%  vs  miR-152  transfected  cells 
29.3%;  MDA-PCa-2b  NC  33.15  vs  miR-152  transfected 
cells  37.59%)  (Figure  6  d,  e,  f).  In  addition  to  reduced 
cell  proliferation,  miR-152  transfected  LNCaP,  PC-3, 
and  MDA-PCa-2b  cells  demonstrated  a  decrease  in  cell 


migration  (Figure  6g)  and  invasion  (Figure  6h). 

Previously,  miR-152  was  demonstrated  to  target 
DNMTl  expression  in  endometrial  tumors  [23].  Since 
we  observed  that  numerous  miRNAs  in  our  panel  were 
influenced  by  methylation,  and  the  methyltransferase 
DNMTl  is  a  regulator  of  DNA  hypermethylation  in 
various  tumor  types,  we  sought  to  explore  the  miR- 
152/DNMTl  relationship  in  PCas.  First,  the  relative 
expressions  of  DNMTl  and  miR-152  were  examined  in 
our  panel  of  cell  lines.  As  expected,  qRT-PCR  showed  that 
DNMTl  expression  was  elevated  in  the  more  aggressive 
cell  lines,  and  this  correlated  with  decreased  expression  of 
miR-152  (Figure7a). 

DNMTl  expression  was  then  examined  in  LNCaP, 
PC-3,  MDA-PCa-2b  cells  transfected  with  miR-152.  Both 
types  of  cells  showed  statistically  significant  decreases 
in  DNMTl  expression  at  the  RNA  and  protein  levels 
(Figure  7b).  This  relationship  is  possibly  direct,  as  miR- 
152  has  putative  binding  sites  in  position  45-54  of  the 
DNMTl  3’  UTR.  To  confirm  this,  we  utilized  a  previously 
established  DNMTl  3'-UTR  luciferase  reporter  system 
[24]  containing  the  miR- 152-binding  sites  (DNMTl  wild- 
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Fig.  4:  Low  miR-152  expression  correlation  with  PCa  metastasis.  Taylor  et  al  data  set  samples  were  sorted  by  average  miR-152 
expression  and  the  mean  from  two  probes  was  normalized  with  z-score.  (a-d)  miR-152  expression  correlates  with  metastases,  invasion 
events,  and  recurrence- free  survival  in  138  PCa  samples  (97  primary  tumors  and  13  metastases).  (a)  miR-152  levels  decrease  monotonically 
fromnomial,  primary  tumor,  and  metastasis,  P<0. 000 1.  (b)  Intensity  in  heatmap,  with  red  corresponding  to  high  and  blue  to  low  expression, 
with  common  event  that  occur  in  PCa  metastasis,  (c)  Number  of  metastatic  events  separated  by  high,  intermediate,  and  low  miR-152 
expression,  (d)  Low  miR-152  expression  correlated  with  decrease  in  probability  of  biochemical  recurrence  free  survival  p=0.0004  (log- 
rank  test). 


African  American  Caucasian  Patients 


Fig  5:  iniR-152  expression  in  AA  and  CA  matched  normal  tnmor  cohort,  (a)  qRT-PCR  of  miR-152  expression  in  nomal- 
tumor  paired  prostate  tissue  samples.  All  samples  were  normalized  to  the  adjacent  normal  samples.  All  values  under  value  one  represent 
down-regulation  of  miR-152.  (b)  Plot  demonstrating  the  individual  miR-152  expression  values  in  AAand  CA  patients:  Note:  50  %  of  the 
individual  AA  patients  had  significant  lower  miR-152  expression  compared  to  only  35  %  of  CA  patients,  p- values  (Benjamini-Hochberg 
corrected  for  multiple  tests) . 
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type  3'-  UTR)  or  mutating  these  sites  (DNMTl  Mu  3'- 
UTR)  that  contain  the  putative  miR-152  binding  sites. 
All  miR-152  transfected  cells  showed  decreased  DNMTl 


wild-type  3'-UTR  luciferase  activity,  whereas  DNMTl 
Mu  3'-UTR  luciferase  activity  was  not  affected  (Figure 
7c).  Since  we  demonstrated  that  loss  of  miR-152  occurs 


LNCaP 


Days  After  Transfection 


b. 


PC-3 


Davs  After  Transfection 


c.  MDA-PCa-2b 


Days  After  Transfection 


d.  NC 


miR-152 


S 
§ 

LNCaP 

"i 

1M.7M  IM0JK»  3.SOOJDOO 

FL2-A 


4M.073  1JXIO.MIO  2M0XI002.KI0J(I00 

FL2-A 


f. 


MDA-PCa-2b 


2«7.IM  IJWOilCO  2.420.702 

FL2-A 


Boyden  Chamber  Migration 

LNCaP  EC;!  MDA-PCa-2b 


NC  miR-152  NC  miR-152  NC  rniR-152 


Fig  6:  Restoring  miR-152  expression  decreases  cell  proliferation,  migration,  and  invasion,  (a)  Proliferation  of  LNCaP,  (b) 
PC-3,  (c)  and  MDA-PCa2b  cells  that  were  transfected  with  either  30  nM  of  miR-152  mimic  or  mlR-NC  (negative  control)  was  measured 
by  MTT  for  6  days.  Results  shown  are  representative  of  three  independent  experiments  averaged  and  normalized  to  Day  1  ±  s.e.  *p<0.05. 
(d)  LNCaP,  (e)  PC-3,  (f)  and  MDA-PCa-2h  cells  were  analyzed  by  flow  cytometry  after  3  days  of  miR-152  mimic  transfection  (30  nM). 
miR-152  caused  G2-M  arrest  in  cell  cycle  progression  in  PCacell  lines.  Results  shown  are  representative  of  three  independent  experiments, 
(g)  Relative  cell  migration  of  LNCaP  PC-3,  and  MDA-PCa-2b  cells  was  measured  utilizing  Boyden  migration  chambers.  Results  shown 
are  representative  of  two  individual  experiments  perfomied  in  triplicate,  (h)  LNCaP,  PC-3,  and  MDA-PCa-2b  cells  transfected  with  miR- 
1 52  mimic  showed  a  decrease  in  the  number  of  cells  invading  through  a  layer  of  Matrigel  relative  to  cells  treated  with  mlR-NC  (negative 
control).  All  data  presented  are  the  means  of  three  independent  experiments  ±  s.e.  *P<0.05. 
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through  promoter  methylation,  we  asked  if  there  was  a 
reciprocal  relationship  between  miR-152  and  DNMTl. 
Treatment  of  LNCaP,  PC-3,  and  ]VIDA-PCa-2b  cells  with 
DNMT 1  siRNA  caused  statistically  significant  increases  in 
miR-152  expression  (Figure  7d). 

To  examine  the  possible  broader  influence  that  loss 
of  miR-152  has  in  PCas,  we  queried  the  TargetScan  in 
silica  database  to  determine  additional  gene  targets  that 
could  be  regulated  by  miR-152.  Of  the  top  genes,  Rictor, 
TGF-P,  SOSl,  ABCD3,  SMAD4,  SOX2,  E2F1,  and  Dicer 
were  predicted.  To  determine  their  influence,  a  custom 
gene  array  of  these  genes  was  designed  (Supplemental 
Table  5),  and  expression  levels  in  LNCaP,  PC-3,  and 
MDA-PCa-2b  cells  that  were  transfected  with  miR-152 
were  assayed.  Although,  the  effect  of  ectopic  miR-152  on 
gene  expression  varied  among  the  cell  lines,  we  observed 
consistent  decrease  expression  in  Rictor,  TGF-P,  and 


SOSl.  Interestingly,  in  the  African  American  derived 
MDA-PCa-2b,  we  observed  significant  decreases  in 
ABCD3  and  SOSl,  which  have  been  previously  reported 
to  be  associated  with  African  American  prostate  tumors  [4, 
5].  Together,  these  results  indicate  that  miR-152  regulates 
epigenetic  events  that  promote  tumorigenesis. 

DISCUSSION 

Numerous  studies  have  now  reported  gene 
differences  AA  and  CA  prostate  tumors  [4,  6,  7,  25,  26]. 
As  the  held  of  health  disparities  is  in  its  infancy,  these 
reports  provide  the  first  evidence  of  population-based 
genetic  contributions  to  aggressive  disease.  However, 
the  regulation  of  these  genes  is  still  elusive.  To  address 
this  question,  without  influence  of  non-epithelial 
stromal  cells  typically  present  in  tissues,  we  utilized 
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Fig.  7:  Comparison  of  miR-152  and  DNMT-1  relationship  in  PCa  cell  lines  (a)  Relative  miRNA  and  mRNA  expressions  of  miR-152 
and  DNMTl  were  determined  by  qRT-PCR  in  a  panel  of  prostate  cell  lines  with  increasing  aggressiveness,  (b)  miR-152  or  scrambled 
oligonucleotides  as  NC  (negative  control)  treatment  of  LNCaP,  PC-3,  MDA-PCa-2b  cells  resulted  in  decreases  of  DNMTl  at  both  RNAand 
protein  levels,  (c)  Dual-luciferase  assays  were  perfonned  for  LNCaP,  PC-3,  and  MDA-PCa-2b  cells  co-transfected  with  the  firefly  luciferase 
constructs  containing  the  DNMTl  wild-type  or  Mu  3'-UTR  and  miR-152  mimics  or  (NC)  negative  control,  (d)  siRNA-DNMTl  treatment 
resulted  in  increases  in  miR-152  expression  in  LNCaP,  PC-3,  and  MDA-PCa-2b  cells,  as  deteraiined  by  qRT-PCR.  DNMTl  expression  was 
analyzed  by  immunoblots  utilizing  anti-DNMTl  antibody  for  both  LNCaP,  PC-3,  a  nd  MDA-PCa-2b  cells  after  siR-DNMTl  treatment. 
For  all  qRT-PCR  experiments,  expression  was  normalized  to  RNU48  (miRNA)  and  GAPDH  (mRNA)  controls.  All  data  presented  are  the 
means  of  three  independent  experiments  ±  s.e.  *P<0.05. 
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AA-  and  CA-derived  cell  lines  derived  from  benign  and 
primary  tumors,  along  with  well-established  cell  lines 
derived  from  metastatic  PCas.  Our  findings  derived  by 
miRNA  microarray  analysis  revealed  a  distinct  miRNA 
expression  pattern  in  AA-derived  cell  lines.  qRT-PCR 
validation  of  these  miRNAs  showed  that  the  five  most 
significant  miRNAs  were  decreased  in  metastatic 
cell  lines,  suggesting  a  role  for  these  molecules  in 
tumorigenesis.  Since  previous  reports  have  demonstrated 
inheritable  epigenetic  changes  in  AA  patients,  we  sought 
to  determine  if  decreased  expression  was  the  result  of 
DNA  methylation.  In  silica  analysis  revealed  that  these 
miRNAs  have  dense  CG  regions  in  the  promoter.  Further 
evidence  that  methylation  of  CpG  islands  causes  silencing 
is  that  each  of  the  miRNAs,  albeit  to  varying  degrees,  is 
re-expressed  after  treatment  with  5-aza-2’d  or  TSA.  To 
our  knowledge,  this  is  the  first  report  that  demonstrates  a 
methylated  miRNA  profile  that  can  be  associated  with  the 
progression  of  AA  prostate  tumors. 

Of  the  miRNAs  influenced  by  5-aza-2’d,  several 
have  been  found  in  other  tumor  types,  including  PCa. 
For  example,  miR-132,  which  targets  both  heparin¬ 
binding  epidermal  growth  factor  and  TALIN2,  is  silenced 
by  hypermethylation  [15].  Of  the  miRNAs  we  assayed, 
miR-152  demonstrated  the  most  significant  reversal  of 
expression,  an  effect  consistent  for  both  LNCaP  and  PC-3 
cell  lines.  Previous  reports  regarding  endometrial  cancer 
[23],  gastrointestinal  cancer  [27],  and  ovarian  cancer  [28, 
29],  and  a  report  relating  to  PCa  [30],  published  during 
the  preparation  of  this  manuscript,  show  that  miR-152 
levels  are  decreased  in  more  advanced  tumors.  Analysis 
of  patient  data  from  the  Taylor  et  al.  study,  deposited  in 
the  NIFI-sponsored  Geo  database,  confinns  that  miR-152 
expression  is  low  in  the  more  advanced  primary  tumors, 
with  the  lowest  expression  in  metastatic  samples.  However 
our  analysis  provides  evidence  that  low  miR-152  levels 
decreases  probability  of  biochemical  recurrence  free 
survival  in  patients.  qRT-PCR  to  measure  normal-tumor 
expression  ratios  confirmed  these  findings,  with  67%  of 
patients  displaying  lower  miR-152  expression.  However, 
the  most  relevant  observation  is  that  AA  patients  display 
decreased  expression  of  miR-152  in  both  uninvolved  and 
paired  tumors  compared  to  CA  patients  of  similar  age, 
stage,  and  Gleason  grade.  To  determine  the  mechanism 
for  decreased  miR-152,  sodium  bisulfite  modification  and 
sequencing  of  miR-152  promoter  in  aggressive  PCa  cell 
lines  were  performed.  Both  LNCaP  and  PC-3  cells  showed 
100%  methylation  status,  highlighting  that,  in  aggressive 
tumors,  hypermethylation  of  the  miR-152  promoter  is  a 
mode  of  silencing  miR-152. 

miR-152  directly  targets  the  3'  UTR  of  DNMTl 
[23].  Across  our  cell  line  panel,  DNMTl  and  miR-152 
showed  an  inverse  relationship  in  expression.  Although 
the  miR-152/DNMTl  inverse  expression  relationship 
was  significant  in  MDA-PCa-2b  cells,  both  the  LNCaP 
and  PC-3  cells  showed  much  greater  miR-152/DNMTl 


expression  differences.  However  in  all  three  cell  lines, 
we  confirmed  that  forced  expression  of  miR-152  results 
in  decreased  expression  of  DNMTl  at  both  the  RNA 
and  protein  levels.  We  also  confirmed  in  LNCaP,  PC-3, 
and  MDA-PCa-2b  cells  decreased  expression  of  several 
candidate  targets  that  have  been  implicated  in  AA  prostate 
tumors.  SOS-1,  which  is  a  regulator  of  EGFR  expression 
and  downstream  signaling,  is  increased  in  AA  PCas,  and 
most  significantly  decreased  in  the  AA  derived  MDA-PCa- 
2b  cells  [5].  Although  Rictor,  a  subunit  of  the  mTORC2 
complex,  has  not  been  directly  implicated  in  AA  PCas,  the 
phosphatidylinositol-3 -kinase/ AKT  (PI3K/AKT)  pathway, 
which  is  a  regulator  of  the  mTORC2  complex,  has  been 
associated  with  AA  PCas  [25].  Whereas  TGF-a  is  a  direct 
target  of  miR-152  in  PCa  cell  lines  [30],  we  observed 
decreased  TGF-P  mRNA  expression  after  forced  miR-152 
expression  in  all  three  cell  lines,  which  further  suggests 
that  loss  of  miR-152  promotes  increased  aggressiveness 
through  multiple  signaling  pathways.  Previously,  we  have 
demonstrated  that  ABCD3  expression  is  associated  with 
African  American  Prostate  tumors  [4.]  Interestingly,  we 
observed  decreased  ABCD3  expression  in  the  MDA-PCa- 
2b  cells  after  miR-152  ectopic  expression,  that  was  not 
observed  in  LNCaP  or  PC-3  cells.  Because  DNMTl  has 
been  characterized  as  one  of  the  main  enzymes  responsible 
for  maintenance  of  global  methylation  patterns  in  tumor- 
related  genes  [31,  32],  we  investigated  the  role  of  DNMTl 
on  miR-152.  We  found  that  depletion  of  DNMTl  through 
siRNA  resulted  in  increased  miR-152  in  both  LNCaP, 
PC-3,  and  MDA-PCa-2b  cell  lines,  highlighting  that  the 
effect  is  independent  of  androgen  sensitivity.  However, 
there  appears  to  be  a  feed-forward  loop  where  either 
loss  of  miR-152  and/or  increased  DNMTl  maintains 
methylation  patterns,  and  hence  miR-152  expression.  miR- 
152  is  gaining  interest  as  a  factor  in  various  tumor  types; 
a  recent  report  concerning  nickel  sulfide-transformed 
human  bronchial  epithelial  (16HBE)  cells  demonstrated 
that  the  miRNA- 152/DNMTl  relationship  develops  early 
in  transformed  cells  [24].  Treatment  of  cells  with  5-aza- 
2’-d  or  depletion  of  DNMTl  led  to  increased  miR-152 
expression  by  reversal  of  promoter  hypermethylation. 
These  findings  are  similar  to  those  presented  here, 
however  the  available  data  relating  to  PCas  suggest  that 
the  relationship  between  loss  of  miR-152  and  increased 
DNMT 1  occurs  during  the  progression  to  advanced  tumor 
status. 

Our  finding  of  a  regulatory  mechanism  that 
maintains  methylation  patterns  has  implications  for  AA 
PCa  patients.  The  fact  that  miR-152  expression  is  lower  in 
non-malignant  tissues  from  AA  patients  compared  to  CA 
non-malignant  tissues  suggests  that  miR-152  expression 
is  regulated  by  inheritable  differences.  Although  the  cause 
of  this  decreased  expression  is  still  speculative,  using 
NCBI  Genome  and  dpSNP  Database  (http://www.ncbi. 
nlm.nih.gov/projects/SNP/snp_ref  cgi?rs=12940701), 
we  did  identify  38  single-nucleotide  polymorphisms 
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(SNPs)  on  chromosome  17q21  (46114527 . 46114613, 

complement!)  region  where  miR-152  is  located.  Of  these, 
rs2001 14569,  which  is  associated  with  a  C  to  G  nucleotide 
change,  and  rs  12940701,  which  is  associated  with  a  C  to 
T  change,  were  of  interest  because  they  may  affect  the 
methylation  status  of  the  miR-152  gene.  The  rsl2940701 
SNP  has  a  relatively  high  frequency  (15%)  of  the  minor 
T  allele  in  the  European  American  population  (n=120) 
relative  to  the  Yoruba  (Nigerian)  population  (n=118). 
Thus,  the  European  American  population  could  have 
lower  methylation  rates  since  the  C  nucleotide  is  changed 
to  the  T  nucleotide,  whereas,  in  the  Yoruba  population, 
this  shift  is  less  frequent,  allowing  for  more  C  nucleotides, 
thus  possibly  increasing  the  rate  of  methylation.  Further 
analysis  should  focus  on  the  methylation  frequency  of 
miR-152  and  on  genotyping  SNPs  that  compromise  miR- 
152  expression  in  AAs  and  CAs. 

In  summary,  the  present  work  highlights  miR-152 
as  a  tumor  suppressor  that  is  inactivated  by  methylation. 
Because  a  large  number  of  tumor/metastasis  suppressor 
genes  are  silenced  as  a  result  of  methylation,  miR-152 
could  be  a  central  regulator  of  key  events  that  contribute 
to  tumorigenesis  and  aggressiveness  and  thus  has  the 
potential  to  be  a  therapeutic  agent  for  PCa  treatment. 
These  results  begin  to  unravel  the  molecular  mechanism 
associated  with  the  aggressive  tumors  of  AA  patients. 

MATERIALS  AND  METHODS 


PCa  cell  lines  and  primary  tissue  samples 

Immortalized  PCa  cell  lines  RC-77T/E  (T3c  poorly 
differentiated  primary  tumor),  RC-77N/E  (non-malignant) 
were  derived  from  AA  patient  as  previously  described 
[21].  RC-43T/E  (T4  poorly  differentiated  primary  tumor) 
and  RC-43N/E  (non-malignant)  were  also  obtained  from 
an  AA  patient  (characterization  unpublished).  Benign  cell 
lines  derived  RC-165N/hTERT  (derived  from  AA  patient) 

[33]  and  RC-170N/hTERT  (derived  from  CA  patient) 

[34] ,  and  primary  tumor  cell  line  RC-92a/hTERT  (derived 
from  CA  patient)  cells  [35].  All  cell  lines  were  cultured  in 
keratinocyte  seram-free  medium  (KGM,  LifeTechnologies, 
Carlsbad,  CA)  supplemented  with  bovine  pituitary  extract, 
recombinant  epidermal  growth  factor,  and  1%  penicillin- 
streptomycin-neomycin  were  maintained  in  KGM  medium 
as  previously  described  [36].  Non-malignant  CA  prostate 
epithelial  cells  (PrEC)  were  obtained  from  Clonetics 
Lonza  (Switzerland)  and  maintained  in  Prostate  Epithelial 
Cell  Growth  Medium  (Clonetics).  The  CA  androgen- 
independent  and  metastatic  PC-3  and  DU- 145  PCa  cell 
lines  were  maintained  in  Dulbecco’s  Modified  Eagle 
medium  (DMEM)  supplemented  with  10%  fetal  bovine 
serum  and  1%  penicillin-streptomycin-neomycin.  The  CA 
LNCaP  and  C4-2b  androgen-dependent  and  -independent. 


respectively,  prostate  cells  were  maintained  in  T-medium. 
RWPE-1  normal  prostate  cells  were  maintained  in  KGM 
as  previously  described  [37].  The  malignant,  androgen- 
receptor  positive  AA  MDA-PCa-2b  cells  were  purchased 
from  ATCC  (Manassas,  VA).  These  cells  were  maintained 
in  F-12K  medium  supplemented  with  20%  fetal  bovine 
serum,  25  ng/ml  cholera  toxin,  10  ng/ml  mouse  epidermal 
growth  factor,  5  pM  phosphoethanolamine,  100  pg/ml 
hydrocortisone,  45  nM  selenious  acid,  and  5  pg/ml  bovine 
insulin. 

Patient  Samples 

Patient  samples  were  obtained  from  the  Cooperative 
Fluman  Tissue  Network  at  the  University  of  Alabama  at 
Birmingham  (UAB)  under  Institutional  Review  Board- 
approved  protocols.  Additionally,  the  Institutional  Review 
Boards  of  Tuskegee  University  and  UAB  approved  the 
use  of  tissues  for  this  study.  Tumor  sections  were  macro 
dissected.  There  was  no  significant  difference  between 
the  two  groups  with  respect  to  tumor  content.  (Average 
percentages  of  tumor  for  AA  and  CA  biopsies  were  both 
54%.)  Paired  normal  tissues  were  also  collected  from  the 
same  patients.  Altogether,  samples  from  each  of  the  20  AA 
and  1 9  CA  patients,  along  with  patient-matched  normal 
adjacent  tissues,  were  used  for  qRT-PCR  validation, 
p-values  were  generated  using  graph  pad  software. 
T-tests  (the  Benjamini-FIochberg  corrected  for  multiple 
tests  that  reflects  the  individual  patients  with  significant 
difference  in  mirl52  expression  relative  to  global  mean 
of  u48  for  Ct  values  and  global  mean  for  delta  Ct  values) 
were  perfonned  to  test  for  significant  differences  between 
the  means,  and  F-tests  were  used  to  assess  differences  in 
variance  between  experimental  groups. 

Quantitative  Real-time  PCR  (qRT-PCR) 

Total  RNA  extraction  was  performed  using  the 
Ambion  recover  all  nucleic  acid  isolation  kit  (AM  1975) 
modified  by  replacing  filters  with  (AM10066G).  RNA 
(10  ng)  was  reverse  transcribed  using  TaqMan  miRNA 
reverse  transcription  kits  (Life  Technologies).  For  inRNA 
expression,  1  pg  of  total  RNA  was  reverse  transcribed 
using  Fligh  Capacity  cDNA  kits  (Invitrogen).  Relative 
expression  of  miRNAs  and  mRNA  was  quantified 
wifh  fhe  TaqMan  Universal  PCR  Masfer  Mix,  No 
AmpErase  UNG,  wifh  fhe  7500  Fasf  Real-Time  PCR 
sysfem  (Life  Technologies).  Thermal  cycling  condifions 
included  enzyme  acfivafion  for  10  min  af  95  °C,  40 
cycles  of  95  °C  for  15  s,  and  60°C  for  60  s,  according  fo 
provider’s  protocol.  Reverse  transcription  for  mRNA  was 
accomplished  as  previously  described  [22,  38].  SYBR 
Green  reagents  (Invitrogen)  were  used  for  quantitative 
real-time  PCR.  Thermal  cycling  conditions  for  primer 
sequences  used  for  mRNA  detection  are  included  in 
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the  table  below  .  Analyses  for  miRNA  and  inRNA  were 
performed  in  triplicate.  RNU48  miRNA,  GAPDH,  and 
18S  ribosomal  RNA  were  used  as  endogenous  controls. 

Transfection  of  miRNA  Precursor 

LNCaP,  PC-3,  and  MDA-PCa-2b  PCa  cell  lines 
were  seeded  at  2  x  10^  cells  in  six-well  plates  on  the  day 
before  transfection.  Lipofectamine  2000  reagent  (Life 
Technologies)  was  used  to  transfect  5-15  nmol  of  miR- 
152  or  5-15  nmol  of  miR-Negative  Control  #1  in  Opti- 
MEM  (Life  Technologies).  Cells  were  harvested  for  assays 
three  days  after  recovery  in  fully  supplemented  media,  as 
described  previously  [38]. 

Cell  Proliferation  Assays 

LNCaP,  PC-3,  and  MDA-PCa-2b  cells  (2.5  x  102) 
were  plated  in  96-well  plates  in  DMEM.  After  being 
cultured  for  24  hr,  cells  were  transfected  with  miR- 
152  or  Negative  Control  #1  (Life  Technologies)  at  a 
concentration  of  30  nM.  Cell  viability  was  determined  at 
24,  48,  72,  96,  and  120  days  after  transfection.  All  cells 
were  incubated  with  10  pi  of  3-(4,5-dimethylthiazol- 
2-yl)-2,5-diphenyltetrazolium  bromide  (MTT,  Sigma- 
Aldrich,  St.  Louis,  MO)  solution  (5  mg/ml)  in  PBS  for  4 
hr.  Dimethylsulfoxide  (50  pi)  was  added  for  10  min  after 
aspiration  of  the  MTT  solution.  Plates  were  read  at  560 
nm. 

Cell  Migration  and  Invasion  assay 

Cell  migration  and  invasion  was  determined  using 
the  Boyden  chamber  assay.  Briefly,  20,000  cells  were 
plated  in  the  upper  chamber,  with  or  without  Matrigel, 
containing  serum-free  medium  containing  1%  bovine 
serum  albumin  for  24  hours;  this  then  was  replaced  with  a 
serum-free  medium  for  an  additional  24  hours.  The  lower 
chamber  contained  medium  with  10  ng/mL  of  epidermal 
growth  factor  was  used  as  a  chemo-attractant.  The  number 
of  cells  that  migrated  or  invaded  through  the  matrix  over 
a  48  hour-period  was  determined  by  counting  cells  that 
stained  with  crystal  violet  on  the  bottom  of  the  filter.  All 
experiments  were  perfonned  in  triplicate. 

Sodium  Bisulfite  Modification  and  Methylation 
Analysis 

DNA  was  extracted  with  GenElute  Mammalian 
Genomic  DNA  Miniprep  kits  (Sigma-Aldrich) 
following  the  manufacturer’s  instructions.  Subsequently, 
bisulfite  conversion  of  DNA  (2  pg)  was  performed 
with  EpiTect  Bisulfite  kits  (Qiagen,  Valencia,  CA). 
The  miR-152  sequence  was  determined  by  using 


the  UCSC  Genome  Browser  [39].  Methylation  PCR 
amplification  primers  for  sequencing  were  designed 
using  MethPrimer  [40].  Sequence  from  5’  to  3’  -  Primer 
1:  GGGTTAGGGGGAGTAGTTAATTTAG  and  Primer 
2:  ATAAACTCCAAAAACATACCCATCA.  qRT-PCR 
was  performed  on  bisulfite-converted  DNA  ( 1  pi),  using 
the  primers  described  above  covering  two  CpG  regions 
upstream  from  the  miR-152  promoter  DNA  sequence. 
A  second  round  of  amplification  was  performed  on  the 
first  PCR  amplicons  and  subsequently  characterized 
by  electrophoresis  on  1.2%  agarose  gels.  Sequencing 
of  amplicons  was  performed  at  the  UAB  Heflin  Center, 
Birmingham,  AL.  Bisulfite  sequencing  was  analyzed  by 
an  online  DNA  methylation  platform,  Bisulfite  Sequencing 
Data  and  Presentation  Compilation  (http://biochem.jacobs- 
university.de/BDPC/). 

Immunoblots 

Cells  were  harvested  3  days  after  miR-152 
transfection,  washed  with  PBS,  and  lysed  in  NP-40  lysis 
buffer  [10  mM  Tris-Cl  (pH  7.4),  10  mM  NaCl,  3  mM 
MgClj  0.5%  NP-40  (Nonidet  P-40),  0.15  mM  spermine, 
0.5  mM  spermidine,  and  protease  inhibitor  solution] 
and  centrifuged  at  12,000  g  for  10  min  at  4°C.  Protein 
concentrations  were  determined  with  the  BCA  Protein 
Assay  kit  (Thermo  Scientific,  Rockford,  IL)  according 
to  manufacturer’s  instructions.  Protein  samples  were 
separated  on  7.5%  SDS-PAGE  pre-cast  gels  (Thermo 
Scientific).  Protein  levels  were  detected  by  anti-DNMTl, 
anti-mTOR,  and  anti-Rictor  rabbit  polyclonal  antibodies 
(Cell  Signaling,  Danvers,  MA),  and  analyzed  by 
chemiluminescence. 

Treatment  with  5-aza-2’-d 

Cells  (60%  confluent)  were  treated  with  5  pM  5-aza- 
2’d  for  5  days,  with  fresh  media  supplemented  with  5-aza- 
2’d  every  day.  Trichostatin  A  (TSA,  a  histone  deacetylase 
inhibitor)  was  administered  at  a  concentration  of  100 
nM  on  the  last  day  of  treatment.  Cells  were  harvested 
and  assayed  for  miR-152  expression  by  qRT-PCR  as 
previously  described  [21,  22]. 

Flow  Cytometry  Analysis 

Cells  were  harvested  3  days  after  miR-152 
transfection,  washed  with  cold  PBS,  fixed,  and 
permeabilized  with  70%  cold  ethanol  for  propidium  iodide 
staining.  Cell  cycle  analysis  was  performed  with  a  flow 
cytometer  (Accuri,  Ann  Arbor,  MI). 
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miRNA  microarray  samples 

For  miRNA  expression  profiling,  11  unique  cell 
lines  were  utilized;  specifically,  PrEC,  RC-77N/E,  RC- 
77T/E,  RC-43T/E,  RC-43N/E,  RC-165N/hTERT,  RC- 
92a/hTERT,  RC-170N/hTERT,  PC-3,  MDA-PCa-2b, 
and  C42B.  Microarray  hybridizations  were  performed 
in  triplicate  for  each  line  by  Ausurgen,  Austin  TX). 
In  addition,  primary  tumor  samples  (n=9)  were  also 
processed  for  microarray  hybridization  (RC-30T,  RC-33T, 
RC-136T,  RC-139T,  RC-28T,  RC-78T,  RC-143T,  RC- 
145T  and  RC-25T).  Clinical  characteristics  are  listed  in 
(supplemental  table  4).  The  protocols  for  RNA  extraction 
and  miRNA  microarray  hybridizations  were  performed  by 
an  outside  vendor  (Asuragen,  Austin,  TX). 

Raw  data  processing  from  microarrays 

For  data  processing,  the  bioconductor  packages 
pipelined  through  R  interfaces  DaNTE  [41]  and  CARMA 
[42]  were  used.  Pre-processing  and  normalization  was 
completed  with  the  Affy  package.  Two  distinct  array 
platforms  were  used  during  the  course  of  this  project;  cell 
lines  were  hybridized  to  one  platform  and  primary  samples 
were  hybridized  to  another.  The  platforms  differ  in  probe 
content  and  replication.  Arrays  of  identical  platforms  were 
quantile  normalized  to  minimize  inter-array  noise  [43], 
using  the  gcRMA  method  [44].  Subsequently,  replicates 
were  merged  into  median  values  for  statistical  association 
analyses. 

Statistical  analysis  of  normalized  data 

For  all  analyses,  p-values  <0.05  were  considered 
as  significant.  Initial  differential  expression  was  detected 
using  an  ANOVA  linear  model.  Of  the  miRNA  genes  that 
showed  significant  variance,  we  measured  the  specific 
differential  expression  among  categorical  groupings; 
diagnosis,  pathological  stage,  and  race.  Differential 
expression  associated  with  categories  was  measured 
using  a  moderated  T-test  in  the  limma  package  of  R. 
Cells  lines  representing  tumor  and  uninvolved  from 
the  same  individuals  were  grouped  into  two  categories; 
Malignant  or  Non-malignant.  This  grouping  also  included 
established  laboratory  cell  lines  (PREC,  PC3  and  C42b). 
We  incorporated  multiple  hypothesis  testing  to  correct  for 
false  positive  (type  1  errors)  and  false  negative  (type  2 
errors).  Specifically,  we  implemented  adjusted  p-values 
(Bonferroni)  and  false  discovery  rates  (FDR)  using  the 
Benjamini  and  Hochberg  method  [45].  A  moderate  10% 
FDR  cutoff  was  used  when  applicable.  We  also  used 
hierarchical  clustering  analyses  to  determine  groups 
of  differentially  expressed  miRNAs  between  specific 
variables  and  to  isolate  specific  nodes  of  probes  that 


shared  similar  expression  trends. 

We  used  non-parametric  tests  when  comparing 
across  array  platforms  to  validate  the  significance  of 
specific  factors  in  a  given  variable  category.  Specifically, 
a  Kruskal-Walis  test  was  used  for  all  variables.  In  addition, 
for  variables  with  only  2  levels,  a  Wilcoxon  Rank  Sum  test 
was  utilized  for  secondary  validation.  In  cases  where  strict 
fold  change  values  were  of  importance,  the  fold  change 
between  two  factor  levels  was  calculated  using  the  ratio 
of  mean  expression  between  the  variables  over  differences 
in  expression  between  the  variables.  [log(x/y)  =  log(x)  - 

log(y)]. 

To  determine  the  ability  of  significant  probe  profiles 
to  differentiate  samples  into  relevant  categories,  we  used 
a  partial  least  squares  (PLS)  principle  component  test. 
Specifically,  once  probes  were  found  to  be  significantly 
differentially  expressed,  following  hierarchical  clustering, 
we  identified  suitable  nodes  or  combination  of  probes,  for 
each  category.  We  determined  differentially  expressed 
probes  based  on  diagnosis,  stage  and  as  the  predictors  and 
miRNA  expression  values  as  the  response.  These  tests 
both  determined  whether  the  probes  statistically  associated 
with  a  specific  variable  or  had  predictive  value  across 
independent  primary  samples  not  utilized  in  the  initial 
discovery  ANOVA  and  non-parametric  tests.  If  the  largest 
percent  variation  could  be  explained  by  the  variable 
extracted/tested,  and  the  samples  segregated  according 
to  proper  category  variable  annotation,  the  probe  set  was 
considered  to  have  significant  predictive  value  for  the 
specified  categories.  For  these  analyses,  the  DAnTE  1.2 
program  was  utilized. 
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Abstract 

Unfixed  tissue  specimens  most  frequently  are  stored  for  long  term  research  uses  at  either  —80°  C 
or  in  vapor  phase  liquid  nitrogen  (VPLN).  There  is  little  information  concerning  the  effects  such 
long  term  storage  on  tissue  RNA  or  protein  available  for  extraction.  Aliquots  of  49  specimens  were 
stored  for  5-12  years  at  —80°  C  or  in  VPLN.  Twelve  additional  paired  specimens  were  stored  for 
1  year  imder  identical  conditions.  RNA  was  isolated  from  all  tissues  and  assessed  for  RNA  yield, 
total  RNA  integrity  and  mRNA  integrity.  Protein  stability  was  analyzed  by  surface-enhanced  or 
matrix-assisted  laser  desorption  ionization  time  of  flight  mass  spectrometry  (SELDI-TOF-MS, 
MALDI-TOF-MS)  and  nano-liquid  chromatography  electrospray  ionization  tandem  mass  spec¬ 
trometry  (nLC-ESI-MS/MS).  RNAyield  and  total  RNAintegrity  showed  significantly  better  results 
for  —80°  C  storage  compared  to  VPLN  storage;  the  transcripts  that  were  preferentially  degraded 
during  VPLN  storage  were  these  involved  in  antigen  presentation  and  processing.  No  consistent 
differences  were  found  in  the  SELDI-TOF-MS,  MALDI-TOF-MS  or  nLC-ESI-MS/MS  analyses  of 
specimens  stored  for  more  than  8  years  at  —80°  C  compared  to  those  stored  in  VPLN.  Long  term 
storage  of  human  research  tissues  at  —80°  C  provides  at  least  the  same  quality  of  RNA  and  protein 
as  storage  in  VPLN. 

Key  words:  assays,  human  tissues,  proteomics,  RNA,  storage,  temperatures 


The  current  gold  standard  for  preservation  of 
RNA  and  protein  analytes  in  tissue  specimens  is 
snap-freezing,  with  subsequent  storage  in  either 
mechanical  freezers  at  —80°  C  or  in  the  vapor 
phase  of  liquid  nitrogen  (VPLN).  Although  RNA 
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and  protein  can  be  isolated  from  tissue  specimens 
processed  for  routine  histologic  analysis  (formalin 
fixation  and  paraffin  embedding),  the  processing 
can  cause  alterations  in  RNA  and  proteins;  how¬ 
ever,  there  is  little  difference  in  assay  results  when 
mRNA  is  analyzed  by  real-time,  reverse  tran¬ 
scriptase  quantitative  polymerase  chain  reaction 
technology  (RT-Q-PCR)  (Steg  et  al.  2006,  2007). 
Frozen  samples  of  tissues,  however,  are  preferred 
for  some  basic  and  translational  studies,  especially 
genome-wide  sequencing  experiments.  Long  term 
storage  of  samples  sometimes  is  desirable,  because 
the  research  value  of  human  specimens,  especially 


DOI:  10.3109/10520295.2014.904927 


1 


Biotech  Histochem  Downloaded  from  informahealthcare.com  by  Dr.  G.S.  Nettleton  on  06/25/14 

For  personal  use  only. 


cancer  specimens,  increases  over  time  if  needed  for 
longitudinal  clinical  and  outcome  data. 

It  is  considered  a  theoretical  advantage  to 
store  tissue  below  the  glass  transition  phase  of 
pure  water,  because  aqueous  based  chemical  reac¬ 
tions  are  thought  to  cease  at  the  glass  transition 
temperature  (Tg).  Such  conditions  are  achieved  in 
VPLN  (  —  150°  C)  and  by  some  mechanical  freezers. 
Owing  to  logistics  and  expense,  however,  many 
investigators  and  biorepositories  store  samples  in 
mechanical  freezers  at  temperatures  of  —70  to  —90° 

C.  It  is  noteworthy  that  there  is  controversy  about 
the  exact  value  of  pure  water;  recent  studies 
suggest  that  this  is  165°  K  or  — 108°  C,  although  it 
frequently  is  reported  to  be  136°  K  or  -137°  C  (Gio- 
vambattista  et  al.  2004).  Complicating  this  issue  is 
the  fact  that  cells  are  not  filled  with  "pure"  water, 
so  the  practical  values  for  the  water  in  mam¬ 
malian  tissues  are  unknown.  Nevertheless,  enzy¬ 
matic  reactions,  in  general,  are  thought  to  continue 
at  —80°  C  and  cells  do  not  remain  viable  when 
stored  at  -80°  C. 

To  compare  the  relative  effects  of  the  two  time- 
honored  methods  for  storing  frozen  tissue  on  RNA 
and  protein  integrity,  we  performed  a  series  of  anal¬ 
yses  to  assess  the  content  and  quality  of  the  RNA 
transcripts  and  proteins /peptides  on  a  cohort  of 
matched  human  tissue  specimens  stored  for  a  num¬ 
ber  of  years  at  both  —80°  C  and  in  VPLN. 

Materials  and  methods 

Patients  and  specimens 

Our  study  was  approved  by  the  Nationwide 
Children's  Hospital  Institutional  Review  Board.  The 
requirement  for  written  informed  consent  from  fhe 
parficipanfs  was  waived,  because  fhe  sfudy  used 
de-idenfified  specimens. 

Aliquofs  of  49  paired  human  tissue  specimens 
were  stored  at  —80°  C  in  mechanical  freezers  and 
matched  aliquots  were  stored  in  VPLN.  An  addi¬ 
tional  12  paired  specimens  were  obtained  from  fhe 
member  institutions  of  fhe  Cooperafive  Human 
Tissue  Network  (CHTN)  and  were  stored  similarly 
in  different  t5q)es  of  freezers.  The  specimens  were 
remnant  solid  tissues  from  surgeries  thaf  were 
available  affer  clinical  diagnoses  had  been  made. 
Tumors,  normal  appearing  tissues  adjacent  to 
tumors  and  other  non-neoplastic  tissue  specimens 
were  used  (Table  1).  The  selection  of  specimens 
for  fhe  firsf  cohorf  was  based  on  fhe  availabilify 
of  mafched  frozen  fissue  stored  in  VPLN  and  at 
—80°  C.  The  specimens  obtained  from  fhe  CHTN  were 
procured  prospecfively  for  fhis  analysis.  All  tissues 
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were  reviewed  by  pathologists  to  ensure  that  the 
matched  specimens  from  both  storage  methods 
were  the  same  histologically  (see  Specimen  Process¬ 
ing  section). 

The  tissue  procured  from  the  CHTN  were 
obtained  within  2  h  of  removal  from  the  patient, 
split  into  identical  aliquots,  snap  frozen  or  snap 
frozen  in  OCT,  then  shipped  to  a  central  facility  on 
dry  ice  for  storage  at  either  —80°  C  or  in  VPLN  for 

I  year.  All  freezers  used  for  this  study  were  moni¬ 
tored  electronically  and  manually  and  no  failures  in 
temperature  maintenance  occurred  during  the  stor¬ 
age  of  the  specimens. 

Specimen  processing  for  RNA  isoiation 

Specimens  not  initially  frozen  in  OCT  were 
embedded  by  placing  the  frozen  fissue  in  par- 
fially  frozen  OCT,  covering  with  additional 
OCT,  and  freezing  at  —20°  C  immediately  prior 
to  analysis.  All  OCT  embedded  specimens  were 
equilibrated  at  —20°  C,  histologic  sections  were 
cut  using  a  cryostat,  and  the  sections  were  stained 
with  hematoxylin  and  eosin  (H  &  E)  for  histologic 
assessment.  Ten  additional  10  pm  frozen  sections 
were  obtained  and  placed  on  dry  ice  for  isolation 
of  eifher  mRNA  or  protein. 

For  RNA  isolation,  the  samples  were  incubated 
for  10  min  in  1  ml  Tri  Reagent  (Molecular  Research 
Center,  Inc.,  Cincirmati,  OH)  at  50°  C  and  insoluble 
material  was  pelleted  by  centrifugation  for  10  min 
at  3,000  X  g.  Nine  hundred  microliters  of  super- 
natanf  were  used  for  RNA  isolation  according  to 
the  manufacturer's  recommendations.  Each  of  the 
resulting  RNA  pellets  was  dissolved  in  25  pi  RNAse- 
free  water.  RNA  was  quantified  using  a  Nanodrop 
ND-1000  specfrophofomefer  (Nanodrop  Technolo¬ 
gies,  Wilmingfon,  DE). 

Quantification  of  totai  RNA  integrity 

For  microcapillary  electrophoresis  measurements 
of  total  RNA  integrity,  the  Agilent  2100  Bioana¬ 
lyzer  was  used  in  conjunction  with  the  RNA  6000 
Nano  LabChip  kits  (Agilent,  Waldbrorm,  Germany) 
following  fhe  manufacturer's  recommendations. 
Bioanalyzer  electrophoresis  report  files  were 
analyzed  for  RNA  integrity  number  (RIN)  and 
degradation  factor  as  described  earlier  (Auer 
et  al.  2003).  For  RIN  analysis,  2100  Expert  software 
version  B.02.03.SI307  (Agilent)  was  used.  RIN 
reports  total  RNA  integrity  on  a  scale  of  1  fo  10, 
where  1  represents  complete  degradation  and 
10  is  the  highest  level  of  RNA  infegrify.  A  value  of 
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Table  1.  Tissue  types  used  to  assess  RNA  integrity  after  long  term  storage  at  -80°  C  and  in  VPLN 


Tissue 

Storage  period 
(years) 

Tissue 

Storage  period 
(years) 

thyroid  papillary  carcinoma 

10 

choroid  plexus  carcinoma 

8 

lymphoma 

10 

neuroblastoma 

8 

benign  neural  tumor 

10 

colon-non  neoplastic 

8 

dysgerminoma 

10 

lipoma 

8 

ganglioglioneuroma 

10 

hepatoblastoma 

8 

pilocytic  astrocytoma 

10 

neuroblastoma 

8 

PNET 

10 

cellular  fibroadenoma  breast 

8 

embryonal  rhabdomyosarcoma 

10 

neuroblastoma 

8 

ganglioglioneuroma 

10 

skull  mass 

8 

spleen,  normal 

10 

myxoid  neoplasm-malignant 

8 

endodermal  sinus  tumor 

10 

Langerhans  histiocytosis 

8 

pancreatic  cystic/solid  papillary  neoplasm 

9 

neurofibroma 

8 

pilocytic  astrocytoma 

9 

neuroblastoma 

8 

hepatoblastoma 

9 

paraganglioma 

7 

ganglioglioneuroma 

9 

ependymoma 

7 

neuroblastoma 

9 

kidney-  non  neoplastic 

7 

ganglioneuroma 

9 

Burkitt’s  lymphoma 

7 

lymph  node,  follicular  hyperplasia 

9 

meningioma 

6 

spleen,  non-neoplastic 

9 

pilocytic  astrocytoma 

6 

MPNST 

9 

ependymoma 

6 

alveolar  rhabdomyosarcoma 

9 

papillary  thyroid  carcinoma 

6 

ganglioglioneuroma 

9 

Burkitt’s  lymphoma 

6 

infantile  myofibromatosis 

9 

osteosarcoma 

6 

neuroblastoma 

9 

ganglioglioma 

5 

pilocytic  astrocytoma 

8 

1  (lowest  integrity)  was  used  for  samples  where 
no  RIN  values  were  available.  For  degradation 
factor  analysis,  results  of  individual  samples  were 
exporfed  from  fhe  2100  Experf  soffware  as  CSV 
files  and  analyzed  by  Degradomefer  soffware  ver¬ 
sion  1.4.2  (Auer  ef  al.  2003).  Degradafion  factor 
reporfs  fofal  RNA  infegrify  on  a  scale  of  1  fo  100, 
where  1  is  fhe  highesf  level  of  RNA  infegrify  and 
100  is  complete  degradafion.  For  samples  where 
no  degradafion  facfors  were  available,  100  (lowesf 
infegrify)  was  used  for  furfher  calculafions. 

Microarray  expression  profiiing 

Samples  were  used  for  expression  profiling  if 
bofh  of  fhe  following  criferia  were  mef:  1)  af  leasf 
one  of  fhe  paired  samples  showed  RIN  >  1,  and  2) 
af  leasf  one  of  fhe  paired  samples  showed  a  deg¬ 
radafion  factor  <  100.  From  37  paired  samples  fhaf 
mef  fhese  criferia,  25  ng  of  fofal  RNA  per  sample 
were  processed  using  isofhermal  SPIA  Biofin 
Sysfem  (NuGEN  Technologies,  Inc.,  San  Carlos, 
CA)  amplificafion;  2.2  pg  cDNA  resulfing  from 
fhe  amplificafion  were  used  for  microarray  hybrid- 
izafion  (Affymefrix  Human  U133A2.0  GeneChips, 
Sanfa  Clara,  CA).The  U133A2.0  microarray  con- 
fains  more  fhan  22,000  probe  sefs  fo  analyze  18,400 


franscripfs  including  14,500  well-characferized 
genes.  Affer  16  h  hybridizafion  af  45°  C,  washing 
and  sfaining  of  microarrays  was  performed  using  a 
Fluidics  Sfafion  450  (Affymefrix);  GeneChips  were 
scarmed  in  a  GeneChip  Scarmer  3000  (Affymefrix). 
All  steps  of  sample  and  microarray  processing 
were  performed  according  fo  manufacfurer's  rec- 
ommendafions.  CEL  files  were  generated  from  DAT 
files  using  GCOS  software  (Affymetrix).  Probe  set 
expression  estimates  were  calculated  using  RNA 
algorithm  in  Array  Assist  Lite  software  version  3.4 
(Stratagene,  Santa  Clara,  CA). 

Characterization  of  mRNA  integrity 

To  evaluate  mRNA  integrity  of  specimens  stored 
imder  differenf  conditions,  rafios  of  signal  intensi¬ 
fies  were  calculafed  for  probe  sefs  fhaf  measure  fhe 
3'  and  5'  ends  of  two  genes  (GAPDH  and  ACTB). 
Because  firsf  sfrand  cDNA  S5mfhesis  uses  oligo(dT), 
RNA  degradafion  can  cause  reduced  signal  intensi¬ 
fies  af  fhe  5'  end  of  franscripfs  and  fherefore  can 
increase  3':5'  rafios.  To  evaluate  whefher  cerfain 
franscripfs  from  specific  genes  showed  consisfenf 
alferafions  due  fo  storage  condifions,  significance 
analysis  of  microarrays  (SAM)  was  used  (Tusher 
ef  al.  2001). 
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Protein  analysis  by  mass  spectrometry 

Frozen  sections  from  matched  pairs  from  eight 
cases  were  cut  at  5  pm  and  mounted  on  glass  micro¬ 
scope  slides  for  histopathological  analysis.  These 
samples  were  selected  from  various  tumor  speci¬ 
mens  in  which  aliquots  of  the  same  specimen  were 
stored  for  >9  years  at  either  —80°  C  or  in  VPLN. 
The  frozen  sections  of  both  aliquots  of  a  pair  were 
reviewed  by  a  pathologist  to  verify  that  each  ali¬ 
quot  was  equivalent  by  microscopic  examination, 
i.e.,  those  that  contained  approximately  the  same 
content  and  histology  of  the  tumor.  Ten  microm¬ 
eter  sections  then  were  cut  from  the  equivalent 
paired  samples  and  mounted  for  analysis  by  mass 
spectrometry.  The  analysis  was  blinded  concern¬ 
ing  which  aliquot  of  each  specimen  was  stored  at 
one  or  the  other  temperature.  When  the  5  pm  sec¬ 
tions  of  both  aliquots  of  a  pair  were  equivalent  by 
microscopic  examination,  e.g.,  contained  approxi¬ 
mately  the  same  content  of  tumor,  up  to  eight 
10  pm  frozen  sections,  depending  upon  the  tumor 
in  the  respective  frozen  section,  were  scraped  from 
the  slide,  avoiding  the  OCT,  and  were  lysed  with 
100  pi  of  lysis  buffer  (20  mm  Hepes,  pH  8.0,  1% 
Tween  20).  The  supernatants  were  separated  using 
a  high  speed  microfuge  at  10,000  x  g  for  10  min 
Eight  pairs  of  samples  of  the  specimens  that  had 
been  stored  years  were  analyzed  (Table  2). 

The  protein  in  each  lysate  supernatant  was  mea¬ 
sured  using  the  Pierce  BCA  protein  assay  and  an 
equivalent  amount  of  protein  (10  pg)  was  loaded  in 
triplicate  with  each  of  the  aliquots  loaded  randomly 
on  one  of  the  eight  sampling  spots  of  an  IMAC-3 
copper  activated  metal  chip  (Bio-Rad,  Hercules, 


CA).  Each  spot  on  the  chip  was  analyzed  by  sur¬ 
faced  enhanced  laser  absorption/ desorption  time  of 
flight  mass  spectrometry  (SELDI-TOE-MS,  Protein 
Biosystem  II;  Bio-Rad).  In  general,  the  approach 
described  earlier  was  used  for  copper  activation 
and  sample  loading  (Adam  et  al.  2002,  Semmes  et  al. 
2005,  McLerran  et  al.  2008a,b,  Grizzle  et  al.  2005a,b). 
Specifically,  the  blinded  samples  were  processed 
using  a  Biomek  2000  (Beckman  Coulter,  Brea,  CA) 
robotic  sample  preparation  platform  that  diluted 
the  samples  and  loaded  them  on  the  IMAC-3  cop¬ 
per  activated  surface.  The  robotic  system  also  spot¬ 
ted  sinapinic  acid  matrix  on  each  sample.  A  control 
sample  was  loaded  and  analyzed  in  at  least  one  of 
the  eight  wells  of  each  metal  chip;  locations  of  cases 
and  controls  on  chips  were  chosen  randomly  to 
minimize  bias  (Adam  et  al.  2002,  Semmes  et  al.  2005, 
McLerran  et  al.  2008a,b,  Grizzle  et  al.  2005a,b). 

The  metal  IMAC-3  chips  also  were  read  on  a 
MALDI-TOP-MS  (Ultraflex  III,  Briiker  Daltronics, 
Billerica,  MA)  using  an  adapter  plate  made  spe¬ 
cifically  for  the  Briiker  instrument.  All  samples 
were  prepared  according  to  the  manufacturer's 
instructions  (see  above).  Our  method  of  MALDI- 
TOE-MS  analysis  has  been  reported  previously 
(Kojima  et  al.  2008).  Specifically,  1600  shots  were 
acquired  automatically  from  2000-100,000  m/z 
with  high  filtering  set  to  on,  deflector  set  to 
1000  m/z,  and  detector  gain  set  at  1769  V,  using 
smart  beam  technology  at  a  single  (empirically 
determined)  laser  energy,  while  using  a  random 
walk  command.  Euzzy  acquisition  was  turned  on 
to  avoid  summing  poor  spectra.  Elexanalysis  then 
was  used  to  baseline  subtract  all  spectra  using  a 


Table  2.  Pairs  of  surgical  specimens  stored  for  at  least  9  years  and  analyzed  by  mass 
spectrometry 


Case  no. 

Diagnosis 

Condition 

Storage  period 

1A 

Langerhans  histiocytosis 

o 

o 

O 

CO 

1 

10  years 

IB 

Langerhans  histiocytosis 

liquid  nitrogen 

2A 

ependymoma 

O 

o 

o 

CO 

1 

9  years 

2B 

ependymoma 

liquid  nitrogen 

3A 

spleen,  non-neoplasia 

O 

o 

o 

CO 

1 

11  years 

3B 

spleen,  non-neoplasia 

liquid  nitrogen 

4A 

lymphoma 

O 

o 

o 

CO 

1 

12  years 

4B 

lymphoma 

liquid  nitrogen 

5A 

hepatoblastoma 

O 

o 

o 

CO 

1 

10  years 

5B 

hepatoblastoma 

liquid  nitrogen 

6A 

lymph  node,  follicular  hyperplasia 

O 

o 

o 

CO 

1 

11  years 

6B 

lymph  node,  follicular  hyperplasia 

liquid  nitrogen 

7A 

MPNST 

O 

o 

o 

CO 

1 

11  years 

7B 

MPNST 

liquid  nitrogen 

8A 

PNET 

O 

o 

o 

CO 

1 

12  years 

8B 

PNET 

liquid  nitrogen 

4 
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top  hat  approach  followed  by  a  batch  text  output 
processing  script  (written  in  house),  thus  allow¬ 
ing  further  processing  to  be  carried  out  in  MatLab 
(Mathworks,  Natick,  MA). 

Specimens  of  lysate  also  were  analyzed  using 
nano-liquid  chromatography  electrospray  ioniza¬ 
tion  tandem  mass  spectrometry  (nLC-ESI-MS/ 
MS)  as  described  earlier  (Wang  et  al.  2010).  For 
these  experiments,  protein  extracts  were  concen¬ 
trated  and  exchanged  using  equal  volumes  of  50 
mM  ammonium  bicarbonafe  fhree  times  using 
3  K  cutoff  spin  fibers  (EMD  Millipore,  Billerica, 
MA)  and  digested  overnight  with  trypsin  gold 
(Promega,  Madison  Wl)  according  to  the  manu¬ 
facturer's  recommendations.  The  resulting  pep¬ 
tides  were  diluted  to  100  ng/pl  in  0.1%  formic 
acid  and  5  pi  were  subjected  to  nLC-ESI-MS/ 
MS  analysis  using  a  ThermoFisher  Scientific,  Inc., 
Walfham,  MA  LTQ-XL  ion  frap  mass  spectrometer 
equipped  with  a  Thermo  MicroAS  autosampler 
and  Thermo  Surveyor  HPLC  pump.  Nanospray 
source,  and  Xcalibur  1.4  instrument  control  (Ther¬ 
moFisher  Scientific,  Inc.).  Proteins  were  searched 
in  species-specific  subsefs  of  fhe  UniRef  database 
(European  Bioinformatics,  Cambridge,  UK).  Tan¬ 
dem  mass  spectrometry  data  were  converted  to 
mzXML  format  using  instrument-specific  con¬ 
version  soffware  (Insfitute  for  Sysfems  Biology, 
Seattle  WA;  Fred  Hutchinson  Cancer  Center)  and 
run  separately  through  SEQUEST  (ThermoFisher 
Scientific,  Inc.)  with  a  "no  enzyme"  setting  so  that 
non-trypsin  cleavage  sites  were  mapped,  and  also 
with  MASCOT  (Matrix  Science  Inc.,  Boston  MA) 
using  a  "strict  trypsin"  setting  so  that  only  trypsin 
cleavage  sites  were  mapped.  ProteoIQ  (BioInquire, 
Athens,  GA)  was  used  to  combine  tandem  mass 
spectrometry  database  search  results  to  determine 
thresholds,  which  identify  as  many  real  proteins 
as  possible  while  encountering  a  minimal  num¬ 
ber  of  false  posifive  protein  identifications.  The 
numbers  of  unique  peptides  were  calculated  per 
sample  based  on  no  enzyme  and  strict  trypsin 
database  searches. 


Statistical  analysis 

To  evaluate  the  significance  of  differences  between 
RNA  and  protein  integrity  from  samples  stored 
under  the  two  storage  conditions,  a  paired  f-test 
was  performed.  Total  RNA  yield,  RIN,  Degradation 
Factor,  and  3'/5'  ratios  of  GAPDH  and  AGTB  were 
analyzed  for  significant  differences.  Transcripts 
with  the  highest  susceptibility  to  storage-induced 
degradation  were  identified  by  paired  analysis  per¬ 
formed  in  SAM  using  a  false  discovery  rate  (FDR) 
<0.12  and  rurming  100  permutations.  The  list  of 
significant  transcripts  expressed  differentially  pro¬ 
vided  by  SAM  was  used  for  EASE  analysis  (Hosack 
et  al.  2003)  to  identify  preferentially  affected  gene 
ontology  categories.  Results  of  EASE  analysis  are 
reported  when  the  EASE  score  is  less  than  0.05  after 
Benjamin!  correction. 

For  the  MALDI-TOF-MS  analysis,  the  MatLab 
toolbox  was  used  to  align  and  peak-select  mass 
spectra,  thus  producing  a  peak  matrix  file  that  then 
was  applied  to  calculate  mean  intensities  and  coef¬ 
ficient  of  variance.  For  the  LC-ESI-MS  data,  Pro¬ 
teoIQ  (NuSep)  was  used  to  incorporate  the  two 
most  common  methods  for  statistical  validation  of 
large  proteome  datasets;  false  discovery  rate  (FDR) 
calculations  combined  with  peptide  and  protein 
probability  approaches  were  used  (Keller  et  al.  2002, 
Nesvizhskii  et  al.  2003,  Weatherly  et  al.  2005).  The 
cutoff  was  selected  at  less  than  1%. 

Results 

RNA 

Microarray  analysis  was  performed  and  mRNA 
integrity  was  assessed  by  3'/5'  ratios  of  GAPDH 
and  AGTB.  Neither  GAPDH  nor  AGTB  showed 
significant  integrity  differences  caused  by  different 
storage  conditions  (Table  3). 

Based  on  the  assumption  that  certain  transcripts 
could  be  especially  susceptible  to  storage-induced 
RNA  degradation,  SAM  was  performed  to  iden¬ 
tify  transcripts  with  significant  differences  between 


Table  3.  RNA  yield  and  integrity  of  specimens  after  iong  term  storage  at  -80°  C  and  in  VPLN 


O 

o 

O 

00 

I 

VPLN 

p  vaiue 

Mean 

Standard  deviation 

Mean 

Standard  deviation 

yieid 

7100  ng 

6800  ng 

4400  ng 

5400  ng 

0.006 

RIN 

5.8 

2.8 

4.3 

3.3 

0.0002 

DF 

39 

37 

56 

42 

0.003 

3':5'  GAPDH 

2.6 

6.2 

2.2 

5.8 

0.9 

3':5'  AGTB 

52 

49 

54 

55 

0.5 
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paired  samples.  All  44  probe  sets  with  the  highest 
significant  differences  between  the  groups  showed 
higher  signals  in  the  group  of  samples  stored  at 
-80°  C  compared  to  those  stored  in  VPLN  (Table  4). 
We  conclude  that  certain  transcripts  are  especially 
susceptible  to  storage-induced  degradation,  which 
causes  signal  loss  in  the  VPLN  stored  samples.  Gene 
Ontology  analysis  of  the  44  transcripts  identified 
by  SAM  showed  significant  over-representation  of 
specific  biological  themes  (Table  5).  This  means  that 


it  is  very  unlikely  that  the  differences  between  the 
two  storage  conditions  could  have  been  caused  by 
random  variation  between  samples. 

Owing  to  the  unexpected  finding  that  samples 
stored  at  —80°  C  showed  greater  RNA  integrity,  we 
investigated  whether  removing  the  samples  from 
storage  boxes  could  be  a  confounding  factor  that 
affects  RNA  integrity.  When  samples  from  a  storage 
box  were  removed,  the  entire  box  was  stored  for 
several  minutes  on  dry  ice  and  this  could  expose 


Table  4.  Probe  sets  with  significantly  lower  signals  in  specimens  stored  in  VPLN  compared  to  -80°  C  according  to  SAM 
analysis 


Probe  set  ID 


Gene  name 


202087_s_at 

20958 1_at 

211911_x_at 

201272_at 

20961 2_s_at 

201649_at 

208729_x_at 

209059_s_at 

211991_s_at 

202675_at 

214864_s_at 

200725_x_at 

210972_x_at 

206559_x_at 

209036_s_at 

211529_x_at 

211799_x_at 

20961 3_s_at 

208870_x_at 

202201  _at 

201231_s_at 

213366_x_at 

202746_at 

200737_at 

210460_s_at 

212085_at 

201931_at 

200820_at 

204599_s_at 

209244_s_at 

217933_s_at 

202474_s_at 

204806_x_at 

218893_at 

200663_at 

217972_at 

218232_at 

202343_x_at 

214836_x_at 

217408_at 

215313_x_at 

210547_x_at 

213160_at 

200752_s_at 


cathepsin  L 

HRAS-like  suppressor  3 

major  histocompatibility  complex,  class  I,  B 

aldo-keto  reductase  family  1 ,  member  B1  (aldose  reductase) 

alcohol  dehydrogenase  IB  (class  I),  beta  polypeptide 

ubiquitin-conjugating  enzyme  E2L  6 

major  histocompatibility  complex,  class  I,  B 

endothelial  differentiation-related  factor  1 

major  histocompatibility  complex,  class  II,  DP  alpha  1 

succinate  dehydrogenase  complex,  subunit  B,  iron  sulfur  (Ip) 

glyoxylate  reductase/hydroxypyruvate  reductase 

ribosomal  protein  L10 

T  cell  receptor  alpha  locus  III  T  cell  receptor  alpha  constant 
eukaryotic  translation  elongation  factor  1  alpha  1 
malate  dehydrogenase  2,  NAD  (mitochondrial) 

HLA-G  histocompatibility  antigen,  class  I,  G 
major  histocompatibility  complex,  class  I,  C 
alcohol  dehydrogenase  IB  (class  I),  beta  polypeptide 

ATP  synthase,  H-i-  transporting,  mitochondrial  FI  complex,  gamma  polypeptide  1 
biliverdin  reductase  B  (flavin  reductase  (NADPH)) 
enolase  1,  (alpha) 

ATP  synthase,  H-i-  transporting,  mitochondrial  FI  complex,  gamma  polypeptide  1 
integral  membrane  protein  2A 
phosphoglycerate  kinase  1 

proteasome  (prosome,  macropain)  26S  subunit,  non-ATPase,  4 
solute  carrier  family  25,  member  6 

electron-transfer-flavoprotein,  alpha  polypeptide  (glutaric  aciduria  II) 

proteasome  (prosome,  macropain)  26S  subunit,  non-ATPase,  8 

mitochondrial  ribosomal  protein  L28 

kinesin  family  member  1C 

leucine  aminopeptidase  3 

host  cell  factor  Cl  (VP  16-accessory  protein) 

major  histocompatibility  complex,  class  I,  F 

Isochorismatase  domain  containing  2 

CD63  antigen  (melanoma  1  antigen) 

coiled-coil-helix-coiled-coil-helix  domain  containing  3 

complement  component  1,  q  subcomponent,  alpha  polypeptide 

cytochrome  c  oxidase  subunit  Vb 

HRV  Fab  N8-VL 

mitochondrial  ribosomal  protein  S18B 
major  histocompatibility  complex,  class  I,  A 
Islet  cell  autoantigen  1 ,  69  kDa 
dedicator  of  cytokinesis  2 
calpain  1 ,  (mu/I)  large  subunit 
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Table  5.  Gene  ontology  categories  over-represented  in  the  44  probe  sets  with  significantiy  lower  signals  in  specimens 
stored  at  VPLN 


System 

Gene  category 

EASE  score 

Benjamin! 

GO  Biological  process 

antigen  presentation 

2X10-® 

3X10-® 

antigen  processing 

2X10® 

3X10-® 

antigen  presentation\,  endogenous  antigen 

4X10-® 

4x  10^® 

antigen  processing\,  endogenous  antigen  via  MHC  class  1 

6X10  ® 

5X10^® 

the  samples  to  a  temporary  temperature  increase. 
The  biorepository's  database  was  used  to  deter¬ 
mine  the  number  of  times  that  a  storage  box  was 
moved  from  VPLN  fo  dry  ice  because  of  roufine 
specimen  refrieval  during  fhe  period  of  sample 
sforage.  Neifher  fhe  number  of  accessions  of  a  box 
nor  fhe  number  of  accessions  of  aliquofs  from  fhe 
same  sample  was  relafed  fo  fhe  loss  of  fhe  fen  mosf 
degradation-sensitive  franscripfs  identified  by  SAM 
(dafa  nof  shown). 

Protein 

The  peaks  generafed  from  fhe  various  aliquofs  ana¬ 
lyzed  by  SELDl-TOF-MS,  which  was  utilized  fo 
assess  whole  "non-enzyme  digesfed"  low  molecu¬ 
lar  weigh!  profeins  and  pepfides,  were  similar  fo 
fhe  equivalenf  friplicafe  measuremenfs  of  each  ali- 
quof  (Fig.  1).  In  some  cases,  fhe  visual  analysis  of  a 
blinded  aliquof  stored  af  —80°  C  yielded  a  higher 
signal  fo  noise  rafio  for  specific  feafures,  while 
fhis  observation  was  reversed  for  ofher  feafures 


depending  on  fhe  f5q)e  of  tissue  specimen  stored 
in  VPLN.  We  could  nof  idenfify  consisfenfly  which 
sforage  condifion  gave  fhe  besf  resulfs  for  analyz¬ 
ing  fhese  samples;  however,  fhe  specfra  obfained 
from  bofh  groups  were  qualifafively  similar  and 
reproducible  despite  fhe  differenf  condifions  and 
fimes  of  sforage. 

Based  on  fhe  qualifafive  resulfs  from  SELDI- 
TOF-MS  analysis,  we  initialed  a  more  exfensive 
analysis  on  fhe  same  samples  using  a  high-end 
MALDI-TOF-MS  insfrumenf  and  a  Briiker  adapfer 
plafe  specifically  made  fo  fif  fhe  SFLDI-MS  probes. 
All  specfra  were  compared  across  each  mafched 
pair  wifh  102  consisfenf  feafures  idenfified  wifh  a 
signal  fo  noise  rafio  >4  (n  =  4).  Similar  resulfs  were 
obfained;  feafures  wifhin  fhe  —80“  C  group  had  a 
coefficienf  of  variation  (CV%)  =0.92,  and  fhe  VPLN 
group  had  a  CV%  =  0.86.  Therefore,  fhe  two  groups 
were  consisfenf  overall.  Specifically,  comparing  fhe 
specfra  of  mafched  pairs,  fhe  majorify  of  peaks  were 
similar  wifhin  each  aliquof  of  fhe  mafched  pairs. 
Some  peaks  were  observed  consisfenfly  af  one 


ecxx) 


8000 


10000 


20  r 


15 


10 


VPLN 


2000 


4000 


6000 


8000 


10000 


Fig.  1.  Spectra  from  a  pair  of  lysed  tissues  from  separate  aliquots  of  the  same  human  specimen  stored  under  different 
conditions  and  for  different  times  at  -80°  C  or  in  VPLN.  The  spectra  are  similar. 
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storage  temperature  (e.g.,  —80°  C),  but  not  at  the 
other  temperature  (e.g.,  VPLN)  and  vice  versa  (i.e., 
absent  at  -80°  C,  but  present  at  VPLN).  Thus,  an 
optimal  temperature  for  storage  could  not  be  deter¬ 
mined  using  either  of  these  approaches  (data  not 
shown). 

Finally,  using  LC-ESI-MS/MS  analyses,  we  com¬ 
pared  the  protein  identification  numbers  (ID's)  of 
matched  pairs  and  a  majority  (183  of  239;  77%)  was 
observed  consistently  in  all  specimens.  The  total 
number  of  proteins  and  the  number  of  protein  ID's 
that  were  independent  of  each  storage  group  also 
were  similar  (total:  independent,  210:27  for  —80°  C 
and  211:29  for  VPLN)  (Fig.  2a).  All  tandem  MS/MS 
spectra  generated  by  the  nLC-ESI-MS  analysis  were 
matched  at  the  peptide  level  against  all  potential 
"no  enzyme"  cleavage  sites  to  assess  differences  that 
likely  are  due  to  either  continued  endogenous  enzy¬ 
matic  degradation  or  shearing  effects  caused  by  the 
storage  method.  This  value  is  estimated  by  assign¬ 
ing  MS /MS  spectra  that  match  peptide  sequences 
with  "strict  trypsin"  cleavage  sites  found  only  at  the 
C-terminal  residue  of  lysine  and  / or  arginine  vs.  "no 
enz5nne"  cleavage  sites  found  at  any  C-  or  N-  terminal 


Comparative  Proteomics 


b)  Comparison  of  Trypsin-strict  To  Random  Peptide  Cleavages 
(VPLN  Vs.  -80°C  Stored  Specimens) 


Fig.  2.  a)  Number  of  unique  peptide  ID’s  generated  by 
LCMS  are  shown  within  a  strict  Venn  diagram  for  the  two 
storage  groups  (-80°  C  and  VPLN)  derived  from  lymphoma 
(L)  specimens  as  an  example,  b)  Number  of  unique 
peptide  identification  numbers  (ID’s)  are  shown  in  the 
figure  for  each  sample  for  all  groups.  L,  lymphoma;  H, 
hepatoblastoma;  F,  follicular  hyperplasia.  Non-neoplastic 
lymph  node  comparing  full  enzyme  (FE)  and  no  enzyme 
(NE)  searches  in  SEQUEST.  In  this  way,  we  can  estimate 
the  number  of  shearing  sites  to  be  expected  from  freezing 
thawing  or  from  a  less  optimal  freezing  approach. 


amino  acid.  While  there  was  a  slight  increase  in 
all  runs  using  the  no  enzyme  vs.  the  strict  trypsin 
approach  (702  vs.  655,  an  approximately  7%  lower 
value  than  702),  there  were  no  consistent  differences 
between  the  two  storage  groups  (Fig.  2b).  Examples 
of  the  most  abundant  proteins  identified  by  mass 
spectrometry  are  listed  in  Table  6. 

Discussion 

Storage  at  —80°  C  and  in  VPLN  are  the  two  most 
common  methods  for  long  term  storage  of  fresh- 
frozen  tissue  samples.  Little  information  is  available 
about  the  relative  benefits  of  these  long  term  stor¬ 
age  methods  for  analysis  of  RNA  or  proteins,  except 
that  storage  near  VPLN  temperatures  is  required  to 
maintain  long  term  viability  of  cells. 

Some  molecules  may  be  affected  by  freezing, 
others  by  thawing  and  others  by  both.  Similarly, 
the  duration  of  long  term  storage  under  either 
condition  probably  results  in  molecular  changes  as 
indicated  by  the  RIN  scores  we  found  for  speci¬ 
mens  whose  average  RINs  were  less  than  expected 
for  specimens  stored  for  a  short  time.  The  freezing 
and  thawing  procedures  for  the  two  storage  meth¬ 
ods  were  equivalent  overall.  The  effects  of  storage 
time  on  RIN  were  based  on  paired  samples  with 
each  tissue  specimen  of  a  pair  stored  for  the  same 
time.  Thus,  the  most  common  difference  between 
each  of  the  paired  samples  and  the  differences  we 
observed  in  RNA  and  proteins  we  attribute  to 
storage  temperature. 

We  found  that  for  49  specimens  of  various  t5q)es 
of  tissues  and  cancers,  storage  at  —80°  C  maintained 
at  least  the  same  RNA  integrity  as  VPLN  storage. 
For  RNA  yield  and  total  RNA  integrity,  —80°  C  stor¬ 
age  provided  significantly  better  results  than  VPLN 
storage.  For  37  specimens  analyzed  by  expression 
profiling,  mRNA  integrity  measured  by  3':5'  ratios 
of  GAPDH  and  ACTB  showed  no  significant  differ¬ 
ences;  by  contrast,  SAM  showed  that  a  number  of 
transcripts  seemed  to  be  especially  susceptible  to 
degradation  in  samples  stored  by  VPLN.  These  tran¬ 
scripts  were  significantly  lower  in  samples  stored  in 
VPLN  compared  to  storage  at  -80°  C.  Transcripts 
that  are  highly  sensitive  to  VPLN  storage  frequently 
are  involved  in  antigen  presentation  and  processing, 
which  suggests  that  these  transcripts  may  be  useful 
for  assessing  storage-induced  RNA  degradation. 

With  the  development  of  real-time  qualitative 
PCR  (RT-Q-PCR),  the  effects  of  degradation  of  RNA 
may  be  minimized  if  short  amplicons  (^  100  bp)  are 
used.  Using  this  technology,  gene  expression  even 
in  paraffin  blocks  can  be  analyzed  reliably  (Steg 
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H,  hepatoblastoma;  LP,  lymphoma;  S,  non-neoplastic  spleen;  rsd,  relative  standard  deviation. 
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et  al.  2006,  2007).  Use  of  RT-Q-PCR  also  reduces 
some  concern  about  long  term  storage  of  tissues 
below  -80°  C. 

We  reported  earlier  that  multiple  freeze-thaw 
cycles  cause  degradation  of  proteins  in  samples  of 
serum.  We  demonstrated  also  that  storage  at  —80° 
C  for  18  months  does  not  cause  significant  degrada¬ 
tion  of  many  proteins;  however,  storage  of  serum 
samples  at  —20°  C  for  more  than  6  months  does 
cause  significant  protein  degradation  (Grizzle  et  al. 
2005a,b).  We  have  shown  here  that  in  the  absence 
of  thawing,  there  were  no  appreciable  differences 
in  protein  integrity  between  the  two  storage  meth¬ 
ods.  While  these  approaches  to  mass  spectrometry 
do  not  evaluate  proteins  and  peptides  that  are  pres¬ 
ent  at  low  concentrations,  for  most  of  the  proteins 
we  evaluated,  there  were  no  differences  between 
the  two  storage  temperatures;  however,  it  should 
be  emphasized  that  the  proteins  analyzed  by  SEL- 
Dl-MS,  MALDI-MS,  and  LC-ESI-MS/MS  represent 
only  a  sample  of  some  of  the  more  abundant  pro¬ 
teins  of  the  proteome.  Our  study  was  limited  to 
global  low  resolution  studies  in  which  the  amounts 
of  protein  were  small  and  the  resources  available 
for  analysis  were  limited;  therefore,  specific  proteins 
that  were  not  observed  may  be  more  or  less  suscep¬ 
tible  to  storage  conditions. 

The  effects  of  sforage  temperature  on  human  tis¬ 
sue  have  not  been  studied  previously  to  the  extent 
that  we  present  here  and  our  results  are  informa¬ 
tive  with  regard  to  the  lack  of  major  changes  in 
the  most  abundant  proteins.  There  is  no  reason  to 
expect  that  changes  during  storage  of  more  abun¬ 
dant  proteins  do  not  mirror  those  of  less  abundant 
proteins.  Together  with  the  RNA  studies,  our  results 
suggest  that  more  extensive,  high  resolution  stud¬ 
ies  are  warranted.  It  is  important  that  most  unique 
proteins  identified  by  mass  spectrometry  in  each 
member  of  the  paired  samples  were  the  same. 
A  relatively  few  unique  proteins  were  not  detected 
in  one  of  the  storage  conditions  compared  to  the 
other  storage  condition  (Eig.  1).  Neither  method  of 
storage  was  clearly  better  than  the  other  with  regard 
to  preserving  specific  proteins  (Table  6). 

Our  results  are  surprising  in  view  of  our  current 
understanding  and  assumptions  regarding  tem¬ 
perature  and  biospecimen  integrity.  Theoretically, 
samples  stored  in  VPLN  should  be  below  the  glass 
transition  temperature  of  water  and  no  enzymatic 
or  hydrolytic  breakdown  of  RNA  or  protein  should 
occur.  We  currently  have  no  hypothesis  about  why 
VPLN  storage  does  not  conserve  RNA  integrity 
as  well  as  —80°  C  storage.  Some  of  the  specimens 
stored  in  VPLN  were  removed  from  the  freezer 
more  frequently  for  distribution  of  specimens  to 


researchers  and  this  could  have  affected  RNA 
integrity  adversely  if  the  specimens  were  warmed 
above  the  glass  transition  phase  for  short  periods 
of  time.  If  this  hypothesis  were  true,  however,  one 
would  expect  that  samples  that  were  moved  more 
frequently  would  show  a  greater  degree  of  RNA 
degradation,  but  we  found  no  such  correlation  in 
our  cohort.  There  is  controversy  concerning  the 
exact  Tg  of  pure  water  (Giovambattista  et  al.  2004) 
and  the  T  of  the  water  phase  of  mammalian  tis¬ 
sues  is  unknown;  however,  enzymatic  reactions  are 
thought  to  continue  at  —80°  C. 

We  conclude  that  the  long  term  storage  of 
fresh-frozen  human  tissue  specimens  in  mechani¬ 
cal  — 80°G  freezers  preserves  at  least  the  same  RNA 
and  protein  quality  as  specimens  stored  in  VPLN. 
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Introduction  &  Objectives:  Prostate  cancer  (PCa)  remains  the  most  frequent  cancer  in  men  in  the 
western  world,  accounting  for  roughly  25%  of  all  new  cases.  Because  of  available  screening  techniques, 
particularly  by  measuring  serum  prostate-specific  antigen  (PSA)  levels,  and  treatment  options,  PCa-specific 
death  represents  only  10%  of  all  cancer-related  deaths.  Treatment  options  are  typically  determined  by  the 
expected  disease  prognosis,  mainly  guided  by  Gleason  scoring  of  a  cancer-positive  biopsy.  Since  only  a 
small  part  of  the  prostate,  typically  less  than  1%,  is  examined  through  histopathology  prior  to  treatment  and 
because  of  fear  for  overtreatment  and  associated  risks,  more  accurate  and  objective  prognosis  is  needed  to 
improve  patient  management. 

Material  &  Methods:  One  hundred  and  two  men  underwent  a  12-core  prostate  biopsy  under  routine  clinical 
practice.  All  patient  data  was  anonymized  and  institutional  review  board  approval  was  obtained.  No  evidence 
of  PCa  was  found  in  20  (20%)  patients.  For  46  (45%)  men,  the  highest  observed  Gleason  score  (GS)  was 
3-1-3,  making  these  potential  candidates  for  active  surveillance  programs.  The  remaining  patients  all  had 
higher-grade  disease,  of  which  26  (25%)  with  GS  3-f4,  8  (8%)  with  GS  4-1-3,  1  (1%)  with  GS  7  and  1  (1%) 
with  GS  4-1-4.  All  individual  biopsy  cores  were  epigenetically  profiled  using  a  quantitative  methylation-specific 
PGR  assay  determining  the  DNA-methylation  status  of  three  markers,  i.e.  GSTPl,  RASSFl  and  APC, 
previously  shown  to  be  associated  with  PCa  development. 

Results:  Unsupervised  hierarchical  clustering  was  performed  identifying  a  group  of  men  whose  biopsies 
were  highly  methylated  for  these  three  markers,  of  which  ~75%  had  higher-grade  disease.  Fligh-grade  PCa 
also  appeared  to  exhibit  a  field  effect,  since  the  correlation  between  increased  methylation  signals  and 
higher  GS  prevailed  even  when  those  cores  with  GS  >  7  were  not  taken  into  account.  A  logistic  regression 
model  including  the  epigenetic  signature,  PSA  and  age  indicated  that  the  overall  average  GSTPl  methylation 
signal  was  the  most  significant  risk  factor.  Patients  were  also  stratified  according  to  the  national 
comprehensive  cancer  network  (NCCN)  risk  categories.  The  central  role  for  GSTPl  was  confirmed  for  the 
NCCN  risk  categories  intermediate  (n  =  27)  and  high  (n  =  l).  When  the  NCCN  low  risk  (n=8)  group  was  also 
included  in  the  aggressive  PCa  category,  RASSFl  was  identified  as  the  best  predictor.  Methylated  signal 
cutoffs  were  determined  for  the  two  most  promising  genes,  GSTPl  and  RASSFl,  showing  that  higher-grade 
tumors  were  associated  with  more  positive  cores. 

Conclusions:  Epigenetic  profiling  of  GSTPl,  RASSFl  and  APC  of  individual  biopsy  cores  could  be  used  to 
evaluate  PCa  aggressiveness  in  an  objective  manner  and  potentially  influence  patient  management.  Long¬ 
term  follow-up  is  required  to  show  that  these  correlations  are  independent  of  Gleason  patterns  or  NCCN  risk 
categories. 


